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1 .  INTRODUCTION 


The  purpose  of  this  report  is  to  provide  a  users  manual  and  make  compara¬ 
tive  runs  to  verify  a  computer  program  called  AERODSN.  This  report  includes 
a  step  by  step  procedure  for  input  and  run  instructions,  a  description  of 
each  output  section  and  comparative  plots  to  verify  normal  force  and  pitching 
moment  coefficients.  The  program  is  a  design  tool  for  preliminary  design 
phase  of  missile  development  programs.  The  current  version  is  an  update  of 
AERODSN  described  in  Reference  1.  Changes  in  the  program  include  referencing 
all  coefficients  to  the  reference  area  (input  by  the  user),  and  changing 
the  input  format  to  free  field  read  for  convenience.  The  major  portion  of 
this  program  is  concerned  with  the  aerodynamics  associated  with  adding  wings 
and  tails  or  both  to  a  body  of  revolution.  A  limited  amount  of  experimental 
body  alone  data  is  included  for  simple  configurations  (ogive  cylinder  or 
cone  cylinder);  however,  body  alone  tabular  data  (CNA.Xcp)  can  be  input.  A 
conical  boattail  segment  can  be  included  if  desired.  Runs  from  AERODSN  were 
compared  with  experimental  data  to  establish  a  range  of  applicability.  Data 
base  configurations  such  as  ADVINT,  HIALFA,  DCATFL  and  COPHED  were  used  for 
comparison. 

II.  PROGRAM  DESCRIPTION 

The  program  was  written  to  be  used  in  preliminary  design  of  typical 
missile  shapes.  Its  primary  function  is  to  calculate  the  aerodynamics  of 
wings  or  tails  or  both  added  to  a  cylindrical-shaped  body  using  a  combina¬ 
tion  of  theoretical  and  imperical  data.  A  limited  amount  of  body  alone 
data  (cone  cylinder  and  ogive  cylinder)  are  included  as  an  option;  however, 
the  body  alone  data  can  be  input  if  desired.  The  roll  orientation  is 
assumed  to  be  a  plus  configuration  or  zero  roll.  Wings  and  tails  must  be 
inline.  Dynamic  stability  derivatives  (damping)  are  calculated  at  zero 
angle  of  attack  and  zero  control  deflection.  Angle  of  attack  data  is  out¬ 
put  corresponding  to  the  values  of  alpha  input.  Table  1  summarizes  the 
applicable  ranges  of  independent  variables.  A  description  of  each  subroutine 
can  be  found  in  Reference  1.  A  listing  of  the  program  is  presented  in 
Appendix  A. 

III.  INPUT 

The  following  tables  describe  each  record  input.  Table  2  is  the 
general  input  format.  Table  3  is  a  sample  input  file  for  a  specific  missile 
configuration. 

IV .  RUN  INSTRUCTIONS 

Two  steps  are  required  to  run  AERODSN  on  the  MICOM  Interdata  8/32 
Computer.  These  steps  are:  (1)  Develop  and  input  file  and,  (2)  Run  the 
program  using  the  CSS  procedure  AERODSN  (already  established  on  the  system). 

A.  STEP  1  -  Referring  to  the  Input  section,  develop  an  input  file 
using  the  edit  mode.  Inputs  are  separated  by  blanks 
or  commas  (free  field  read).  Save  the  file  by  typing 
"SAVE  filename".  Exit  the  edit  mode  by  typing  "END". 
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PHSCXailO  BUM 


wuM-Mrt  num 


TABLE  1. 

MACH  NUMBER 
ANGLE  OF  ATTACK 
ROLL  ANGLE 

CONFIGURATIONS 

Sweep  angle 

FIN  TAPER  RATIO 


APPLICABLE  RANGE  OF  INDEPENDENT  VARIABLES 


0.0  TO  4.0 
0  TO  20  DEG. 

0.0  DEG. 

BODY , BODY  TAIL 

BODY  WING, BODY  WING  TAiL 

NO  FORWARD  SWEPT  WINGS. 


0-1. 


ASPECT  RATIO 


CONSISTANT  WITH  SLENDER  BODY  THEORY 


TABLE  2.  GENERAL  INPUT  FORMAT  AND  DEFINITION  OF  INPUTS 


ALL  GEOMETERIC  VARIABLES  MUST  hAVS  THE  SAME  SYSTEM  OF  UNITS. 
I.E.,  CALIBERS,  INCHES,  ETC.  ALL  ARE  DIVIDED  61  DREF  WITHIN 
THE  PROGRAM.  THE  PROGRAM  WORKS  IN  CALIBERS. 

A  CALIBER  IS  DEFINED  AS  1  REFERENCE  DIAMETER. 

INPUT  BODY  CENTER  OF  PRESSURE  (POSITIVE  AFT  OF  NOSE)  -  (CAL) 
PROGRAM  USES  BRITISH  DATA  SHEETS  FOR  BASIC  WING  LIFT  AND  XCP. 
PROGRAM  USES  SLENDER  BODY  THEORY  (NACA  1307)  FOR  WING-BODY 
AND  WING-TAIL  INTERFERENCE. 


ALL  INPUTS  ARE  FREE  FIELD.  INPUTS  CAN  BE  SEPARATED  BY  A  COMMA, 
ONE  OR  MORE  BLANKS,  OR  BY  A  SLASH  (/). 

ALL  VARIABLES  IN  AN  INPUT  LINE  MUST  HAVE  A  VALUE  TO  READ, 

I.E.,  0.  FREE  FIELD  READ  WILL  NOT  INTERPRET  BUNKS  AS  ZERO'S. 

.  INPUT  1 


TITLE  CARD  (1-80  CHARACTERS) 


INPUT  2 


NWPOS 

s 

1 

INPUT  XWING  AS  STATION 

OF  WING  LEADING  EDGE 

n 

s 

2 

PROGRAM  WILL  MOVE  WING 

1/4  CHORD  TO  WING  HINGE 

LINE 

n 

_ 

3 

PROGRAM  WILL  MOV?  WING 

1/2  CHORD  TO  WING  HINGE 

LINE 

»• 

m 

A 

INPUT  XWING  AS  STATION 

OF  WING  TRAILING  EDGE 

NTPOS 

S 

1 

INPUT  XT AIL  AS  STATION 

OF  TAIL  LEADING  EDGE 

it 

: 

2 

PROGRAM  WILL  MOVE  TAIL 

1/4  CHORD  TO  TAIL  HINGE 

LINE 

n 

s 

3 

PROGRAM  WILL  MOVE  TAIL 

1/2  CHORD  TO  TAIL  HINGE 

LINE 

»» ^ 

m 

4 

INPUT  XTAIL  AS  STATION 

5 

OF  TAIL  TRAILING  EDGE 

TABLE  2 


GENERAL  INPUT  FORMAT  AND  DEFINITION  OF  INPUTS  (Cont'd) 


mcdvt 

It 


NOUT 


NRUN 

NDBASE 


=  1  NO  EFFECT. 

=  2  USE  NACA  REPT.  1253  (MCDEVITT)  FOR  FIN  ALONE  CLA. 
GOOD  FOR  RECTANGULAR  FINS  ONLY. 

GOOD  FOR  MACH  NUMBERS  BELOW  1.1  ONLY. 

=  1  TRIM  OUTPUT  ONLY. 

=  2  ALPHA  OUTPUT  ONLY. 

=  3  TRIM  AND  ALPHA  OUTPUT, 
s  INITIAL  RUN  NUMBER.  -  1  TO  4  DIGITS. 

=  1  NO  EFFECT. 

=  2  OUTPUT  DATA  BASE  FILE. 


INPUT  3 


NOSE 

= 

1  OGIVE  NOSE 

n 

s 

2  CONICAL  NOSE 

XLN 

* 

NOSE  LENGTH 

XL 

- 

TOTAL  LENGTH 

INPUT  4 


RBW  s  RADIUS  OF  BODY  AT  WING 

RBT  c  RADIUS  OF  BODY  AT  TAIL 

DREF  s  REFERENCE  LENGTH 

AREF  s  REFERENCE  AREA 

DREF FT  s  REFERENCE  LENGTH  -  FT. 

AREFFT  s  REFERENCE  AREA  IN  SQ.  FT. 

YIY  3  TRANSVERSE  MOMENT  OF  INERTIA  -  SLUG-FT2. 
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TABLE  2.  GENERAL  INPUT  FORMAT  AND  DEFINITION  OF  INPUTS  (Cont’d). 


.  INPUT  5 

NBODY  =  NUMBER  OF  POINTS  FOR  BODY  ALONE  CNACY  AND  XCPCY 

(IF  NB0DY=0,  OMIT  IMPUTS  6.  PROGRAM  WILL  COMPUTE  SLOPE  DATA.  ) 

.  INPUTS  6 

TMACh  =  TABLE  OF  MACH  NUMBERS  FOR  CNACY  AND  XCPCY 

TCNACY  =  TABLE  OF  BODY  ALONE  NORMAL  FORCE  COEFF  SLOPES  (BASED  ON 
BODY  CROSS  SECTIONAL  REFERENCE  AREA)  -  1 /RAD 
TXCPCY  =  TABLE  OF  BODY  ALONE  CENTER  OF  PRESSURE  -  CAL  FROM  NOSE 

.  INPUT  7 

NBTL=0  PROGRAM  WILL  NOT  ADD  BOATTAIL  CNA  AND  XCP. 

"  =1  PROGRAM  USES  BOATTAL  CARDS  TO  CALCULATE  BOATTAIL  EFECTS 

DBOD  =  DIAMETER  RATIO  (BOATTAIL/CYLINDER) 

XLBOD  =  RATIO  (BOATTAIL  LENGTH/CYLINDER  DIAMETER) 

.  INPUT  8 

NMF  =  NUMBER  OF  POINTS  FOR  XCG  VS  MACH  NUMBER 

. . .  INPUTS  9 

TMF  s  TABLE  OF  MACH  NUMBERS  FOR  XCG 

TXCG  s  TABLE  OF  CENTER  OF  GRAVITY  -  CALIBERS  FROM  NOSE. 

TWEIGH  s  TABLE  OF  WEIGHT  -  LBS. 


INPUT  10 
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r ABLE  2.  GENERAL  INPUT  FORMAT  AND  DEFINITION  OF  INPUTS  (Conrluded). 


NFMA  =  NUMBER  OF  FREE  STREAM  MACH  NUMBERS 

.  INPUTS  11 

TFMA  =  TABLE  OF  FREE  STREAM  MACH  NUMBERS 

TALT  =  ALTITUDE  FOR  EACH  MACH  NO.  -  FT  ABOVE  SEA  LEVEL 

THE  =  REYNOLDS  NUMBER  FOR  EACH  MACH  NO.  -  MILLIONS/FT . 

NOTE:  PROGRAM  MULTIPLIES  TRE  BY  i  MILLION  TO  GET  RE. 

.  INPUTS  12-26 


TBOW 

= 

TABLE 

OF 

WING 

EXPOSED  SEMi -SPANS (LOAD 

0.  IF  NO  wing; 

TCRW 

= 

TABLE 

OF 

WING 

ROOT  ChORDS 

TTRW 

= 

TABLE 

OF 

WING 

TAPER  RATIOS(CT/CR) 

TSWTEW 

= 

TABLE 

OF 

wING 

TRAILING  EDGE  SWEEP  ANGLES(DEG) 

TXWING 

r 

TABLE 

OF 

WING 

STATIONS  (MEASURED  FROM 

NOSE)  (SEE  NWPGS) 

T.iINGW 

- 

TABLE 

OF 

WING 

HINGE  LINES 

TBOT 

= 

TAELE 

OF 

TAIL 

EXPOSED  SEMI -SPANS (LOAD 

0.  IF  NO  TAIL) 

TCRT 

s 

TABLE 

OF 

TAIL 

ROOT  CHORDS 

TTRT 

= 

TABLE 

OF 

TAIL 

TAPER  RATIOS  (CT/CR) 

TSWTET 

= 

TABLE 

OF 

TAIL 

TRAILING  EDGE  SWEEP  ANGLES  (DEG) 

TXTA1L 

= 

TABLE 

OF 

TAIL 

POSITIONS  (MEASURED  FROM 

NOSE) 

THINGT 

- 

TABLE 

OF 

TAIL 

HINGE  LINES 

"  =0.  (PROGRAM  SETS  HINGT=1/4  MEAN  AERODYNAMIC  CHORD) 

TDELT  =  TABLE  OF  TAIL  DEFLECTIONS  -  DEG. 

TDELW  =  TABLE  OF  WING  DEFLECTIONS  -  DEG. 

TALP  —  TABLE  OF  ANGLES  OF  ATTACK  -  DEG. 
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TABLE  3.  SAMPLE  INPUT  FILE  FOR  AER0D3N 


■  SAMPLE  INPUT  FILE:  BODY-WING-TAIL. 

2.  1  i  i  3  1  1 

3.  i  7.470  47.000 

4.  3.  3.  6.0  26.21  .50  .1903  1 

5.  b 

6.  .5  0.197  3.014 

/•  .6  3.059  3. 154 

6.  .95  3.24b  3.317 

9.  1.2  3.151  2.431 

10.  1.5  3.  145  1.976 

11.  1.6  3.00b  1.742 

12.  0  i.O  0.0 

13.  2 

14.  0.0  4.544  i.O 

15.  1.0  4.544  1.0 

16 .  2 

17.  .5  0.0  3.592 

16.  1.8  0.0  12.933 

19.  7.149 

20.  3.051 

21.  1.0 

22.  20.0 

23.  32.32 

24.  0.0 

25.  5.974 

26.  3.051 

27.  i.O 

26.  20.0 

29.  41.226 

30.  42.582 

31 .  0. 

32.  0.0 

33.  -5  -2  0.0  2  4  6  8  10  12  14  16  18  20 
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B.  STEP  2  -  To  run  the  program  typ.  "AERODSN"  and  press  "RETURN". 

Wait  for  instructions  to  appear  and  enter  the  name  of 
the  data  file.  When  the  program  has  finished  running 
an  will  appear  on  the  screen. 

V.  oirrrjx 

The  output  consists  of  five  major  groups: 

.  GEOMETRIC  VARIABLES 

.  FLIGHT  CONDITIONS 

.  ZERO  ANGLE  OF  ATTACK  DATA  (SLOPES) 

.  TRIM  POINT 
.  ANGLE  OF  ATTACK  SWEEP 

Presented  in  Table  A  are  definitions  of  each  output  variable  calculated 
within  each  output  group.  Table  5  illustrates  a  sample  output. 

The  group  GEOMETRIC  VARIABLES  is  divided  into  three  components  which 
are  body,  wing  and  tail.  The  corresponding  geometric  values  are  listed  in 
column  form  with  the  appropriate  symbol  ant;  definition.  Printed  under 
FLIGHT  CONDITIONS  are  Mach  number,  center  of  gravity  location  (calibers 
from  the  nose),  altitude  (ASL)  in  feet  and  Reynolds  number  (per  ft.).  These 
values  are  taken  directly  from  the  input  file.  The  ZERO  ANGLE  OF  ATTACK 
DATA  (SLOPES)  group  is  divided  into  body,  wing,  tail  and  total  configura¬ 
tion  values.  Each  aerodynamic  coefficient  slope  is  evaluated  at  zero  angle 
of  attack  and  is  listed  in  column  form  corresponding  to  the  missile  component. 
Aerodynamic  coefficients  are  calculated  for  individual  angles  of  attack  in 
the  ANGLE  OF  ATTACK  SWEEP  group. 

VI.  PROGRAM  VAT  I PATTON 

Comparisons  are  made  to  establish  the  range  of  validity.  Configura¬ 
tions  taken  from  the  AERODYNAMIC  DATA  BASE  such  as  ADVINT,  HIALFA,  DCATFL , 
and  COPHED  were  used.  Normal  force  and  pitching  moment  coefficients  were 
generated  by  AERODSN  and  plotted  against  wind  tunnel  data  for  comparison. 

Detailed  geometric  description  of  each  missile  configuration  tested 
is  included  with  the  plots.  The  Mach  numbers  range  from  0.2  to  4.0. 

Angle  of  attack  ranges  from  -5  to  25  degrees.  Angles  of  attack  greater 
than  25  degrees  are  included  in  this  report;  however,  user  discretion  is 
advised.  Missile  configurations  include  body,  body-tail,  body-wing  and 
body-wing-tail . 

Figures  1  and  2  describe  the  geometry  of  ADVINT.  Figures  3-26  show 
comparison  between  ADVINT  wind  tunnel  data  and  AERODSN.  Normal  force  and 
pitching  moment  coefficients  are  shown  at  Mach  numbers  of  1.95,  3.01  and 
4.02.  Body,  body-tail,  body-wing  and  body-wing-tail  configurations  are 
compared.  Run  numbers  beginning  with  A  are  ADVINT  wind  tunnel  data  and 
runs  beginning  with  G  are  AERODSN  runs. 
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TABLE  4.  DEFINITION  OF  OUTPUT 


FLIGHT  CONDITIONS 

r'MACri  FREE  STREAM  MACH  NUMBER. 

XGG  LOCATION  OF  CENTER  OF  GRAVITY  (CALIBERS  FhOM  NOSE). 

WEIGHT  WEIGHT  (LBS.). 

ALT  ALTITUDE  ABOVE  SEALEVEL  (FT.). 

RE  REYNOLDS  NUMBER  X  lOEo.  (PER  FT) 

ZERO  ANGLE  OF  ATTACK  DATA  (SLOPES) 

EKUP  UPWASH;  RATIO  OF  LIFT  ON  WINGS  IN  PRESENCE  OF  BODY  TO  LIFT 
OF  WING  ALONE. 

BKCC  CARRYOVER;  RATIO  OF  LIFT  ON  BODY  DUE  TO  WING  TO  LIFT  OF 
WING  ALONE. 

CLA  LIFT  COEFFICIENT  OF  TAIL  OR  WING, BASED  ON  RESPECTIVE  FIN  AREAS. 

(1/RAD) . 

CPF  CENTER  OF  PRESSURE  OF  WING  OR  TAIL  MEASURED  FROM  ROOT 

CHORD  LEADING  EDGE. 

CNA  LIFT  COEFFICIENT  SLOPE  OF  COMPONENTS  BASED  ON  AREF  (1/RAD). 

XCP  COMPONENT  CENTER  OF  PRESSURE  DIVIDED  BY  DREF. 

MEASURED  FROM  NOSE. 

CMACG  PITCHING  MOMENT  COEFFICIENT  ABOUT  CG  BASED  ON  AREF  AND 
DREF. 

CMQ  PITCH  DAMPING  COEFFICIENT  SLOPE  (ABOUT  C.G.)  DUE  TO  PITCH  RATE,  BASED 

ON  AREF  AND  DREF  (1/RAD). 

13 
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TABLE  4.  DEFINITION  OF  OUTPUT  (Cont'd) 

CLP  ROLL  DAMPING  COEFFICIENT  SLOPE  BASED  ON  AREF  AND  DREF. 

(1/RAD). 

CLD  ROLL  MOMENT  COEFFICIENT  SLOPE  (4  FINS  OR  4  WINGS  DEFLECTED), 

BASED  ON  AREF  AND  DREF  (1/RAD). 

CHA  HINGE  MOMENT  COEFFICIENT  SLOPE  DUE  TO  ANGLE  OF  ATTACK, 

BASED  ON  AREF  AND  DREF  ( 1 /RAD) . 

CHD  HINGE  MOMENT  COEFFICIENT  SLOPE  DUE  TO  CONTROL  DEFLECTION, 

(FOR  ALL  MOVEABLE  FINS), BASED  ON  AREF  AND  DREF  (VRAD). 

CND  NORMAL  FORCE  COEFFICIENT  SLOPE  DUE  TO  CONTROL  DEFLECTION, 

BASED  ON  AREF  (i/RAD). 

CMDCG  PITCHING  MOMENT  COEFFICIENT  SLOPE  ABOUT  CG  DUE  TO 
CONTROL  DEFLECTIONS  (i/RAD). 

SIGMA  PITCH  OSCILLATION  NATURAL  FREQUENCE  (CYCLES/SEC). 

TRIM  POINT 

ALPHA  TRIM  ANGLE  OF  ATTACK  (DEG). 

DELW  WING  DEFLECTION  (DEG). 

DELT  TAIL  DEFLECTION  (DEG). 

CN  TOTAL  NORMAL  FOHCE  COEFFICIENT  (NF/U  AREF). 

XCP  LOCATION  OF  CENTER  OF  PRESSURE  (  FROM  NOSE)  XCP/DREF. 

CMCG  TOTAL  PITCHING  MOMENT  COEFFICIENT  ABOUT  C.G.  (PM/Q»AREF»DREF) 
HMW  HINGE  MOMENT  COEFFICIENT  OF  WING  HMW/(Q*AREF*DREF) . 

HMT  HINGE  MOMENT  COEFFICIENT  OF  TAIL  HMT(/Q*AREF*DREF) . 

GLOAD  LATERAL  MANEUVERABILITY  (G’S) 
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TABLE  4.  DEFINITION  OF  OUTPUT  (Concluded) 


ANGLE  OF  ATTACK  SWEEP 

ALPHA  ANGLE  OF  ATTACK  (DEGj . 

CN  TOTAL  NORMAL  FORCE  COEFFICIENT  NF/(Q»AREF). 

XCP  LOCATION  OF  XCP  (  FORM  NOSE)  XCP/DREF 

CMCG  TOTAL  PITCHING  MOMENT  COEFFICIENT  ABOUT  C.G.  PM/ (Q*AREF*DREF) 
CHW  HINGE  MOMENT  COEFFICIENT  OF  WING  HMW/ (Q#AREF*DREF) . 

CHT  HINGE  MOMENT  COEFFICIENT  OF  TAIL  HMT/ (Q»AREF*DREF) . 

CNB  BODY  NORMAL  FORCE  COEFFICIENT  NFB/ (Q*AREF) 

XCPB  LOCATION  OF  BODY  XCP  (  FROM  NOSE, XCPB/DREF) . 

CNW  NORMAL  FORCE  COEFFICIENT  OF  WING  PLUS-BODY  INTERFERENCE 

NFW/(U*AREF) . 

XCPW  LOCATION  OF  WING  XCP  (  FROM  NOSE)  XCPN/DREF . 

CNT  NORMAL  FORCE  COEFFICIENT  OF  TAIL  PLUS  BODY  TAIL  INTERFERENCE 

NFT/ (Q#DREF ) . 

XCPT  LOCATION  OF  TAIL  XC?  (  FROM  NOSE)  XCPT/DREF. 

CNTV  NORMAL  FORCE  COEFFICIENT  ON  TAIL  DUE  TO  WING  VORTICES. 
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Figure  27  describes  the  Copperhead  geometry.  Figures  28-33  show 
normal  force  and  pitching  moment  coefficient  generated  by  AERODSN  compared 
to  that  of  Copperhead  wind  tunnel  data.  The  Mach  numbers  tested  are  0.5, 

0.8  and  0.95  for  the  body-wing-tail  configuration.  Runs  beginning  with  C 
are  Copperhead  wind  tunnel  data  and  runs  beginning  with  T  are  AERODSN  runs. 

Figures  34-36  describe  the  geometry  of  the  missile  shape  used  for 
high  angle  of  attack  comparison.  Nose  1,  afterbody  1  and  tail  11  were  used. 
Alpha  ranges  from  25-45  degrees  at  Mach  numbers  of  2.0,  2.5  and  3.01. 

Runs  beginning  with  "B"  are  HIALFA  wind  tunnel  data  and  runs  beginning  with 
"0"  are  AERODSN  runs.  Figures  37-42  show  normal  force  and  pitching  moment 
coefficient  comparison. 

Figures  43-44  describe  the  geometry  for  the  body  tail  configuration 
used  in  this  comparison.  Figures  46-56  show  normal  force  and  pitching 
moment  comparison  at  Mach  numbers  of  0.2  to  1.0.  Angle  of  attack  ranges 
from  -6.0  to  6.0.  Runs  beginning  with  "A”  are  DCATFL  wind  tunnel  data  and 
runs  beginning  with  "E"  are  AERODSN  runs. 

VII.  SUMMARY  AND  CONCLUSIONS 

The  preliminary  design  computer  program  "AERODSN"  has  been  evaluated 
for  accuracy  and  range  of  applicability.  A  description  of  the  inputs  and 
outputs,  as  well  as  run  instructions  for  running  the  program  on  MICOM's 
Interdata  8/32  computer,  is  included.  "AERODSN"  was  tested  at  Mach  numbers 
from  0.0  to  4.0.  All  configurations  were  tested  with  an  angle  of  attack 
range  from  0.0  to  20.0  degrees  with  the  exception  of  one  tested  from  25.0 
to  45.0  degrees.  Roll  orientations  are  zero  for  all  cases.  Body,  body-wing, 
body-tail  and  body-wing-tail  configurations  were  tested  through  the  alpha 
and  Mach  ranges  indicated  based  on  the  comparative  analysis.  The  following 
conclusions  are  drawn: 

A.  The  applicable  Mach  range  is  from  0  to  4. 

B.  Tn  general,  good  results  were  obtained  in  CN  and  CM  at  angles  of 
attack  up  to  20  deg  alpha.  In  some  cases,  where  configurations  depart 
from  the  missile  shapes  applicable  for  this  program,  comparisons  are  not 
good  for  angles  of  attack  greater  than  5  deg. 

C.  "AERODSN"  is  applicable  for  configurations  with  roll  orientations 
of  zero  only. 

D.  "AERODSN"  is  an  acceptable  preliminary  design  tool  for  obtaining 
a  quick  aerodynamic  profile  for  missile  shapes  described  in  this  report. 
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Figure  15.  Body  wing  normal  force  vs  angle  of  attack  at  M  =  1.05. 
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vs  angle  of  attack  at  M  =  4.02. 
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Figure  21.  Body  wing-tail  norma)  force  coefficient  vs  angle  of  attack  at  M  =  1.93. 
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Figure  31.  Pitching  moment  coefficient  vs  angle  of  attack  at 
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Figure  36.  Details  of  straight  fin  configurations. 
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Figure  40.  Pitching  moment  coefficient  vs  angle 
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Figure  44.  Fin  details. 
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Figure  46.  Pitching  moment  coefficient  vs  angle  of  attack  at 
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Figure  48.  Pitching  moment  coefficient  vs  angle  of  attack 
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Figure  51.  Normal  force  coefficient  vs  angle  of  attack  at 


SYM  RUN  CONFIGURATION 


ANGLE  OF  ATTACK ,  ALPP,  DEG 


ANGLE  OF  ATTACK,  ALPP,  DEG 


RUN  CONFIGURATION 


E 


Figure  55.  Normal  force  coefficient  vs  angle  of  attack 


CONFIGURATION 


Pitching  moment  coefficient  vs  angle  of  attack 


APPENDIX  A 
PROGRAM  LISTING 


SBATCM 

PHELImInAK Y  UfcSiGN  PHUGK AM  iltH  PAKAHfciEKIC  »AKIAfiU». 

ALL  GEUmFTERIC  VARIABLES  MUSI  HAVE  THE  SA«F  SYSTEM  UP  UNITS. 
I.E.,  CALIBERS,  INCHES,  ETC.  ALL  AHt  UivlDF.D  BY  OH t F  mITmIn 
THE  PROGRAM,  XrtE  PROGRAM  «URKS  IN  CALlHEKS. 

A  CALI bEK  IS  DEFINED  AS  1  KEfERENCE  DIAnETEH. 

INPUT  BODY  CENTER  UF  PKESSUHE  (PUSITlVt  AFT  OF  NOSE)  -  CCAL) 
PROGKAM  USES  BRIT ISH  DATA  SHEETS  FUH  BASIC  wlNG  LIU  Axu  aCP. 
PROGkAM  USES  SLENDER  BUUY  THEORY  (NACA  1107)  FUH  «lxG-HJO( 

AND  wlNG-TAlL  INTERFERENCE. 


ALL  Inputs  are  free  FIELD.  INPUTS  CAN  he  SEPAkaIED  BY  a  COMMA, 

One  UK  MURE  BLANKS,  UR  BY  A  SLASH  (/). 

ALL  VARIABLES  IN  AN  INPUT  Ll've.  MUST  HAvE  A  VALUE  lu  hEAU, 

I  .£.  ,  U.  FREE  FIELD  READ  » ILL  NUT  IMEKPHET  BLANKS  AS  EERO'S. 

. . . . . . .  input  t 

IITLE  CARD  (l-«0  CHARACTERS) 

. . . . . . .  INPUT  i 

n»pu5  *  l  input  xwing  as  station  uf  wing  leading  edge 

"  *  2  PROGRAM  «ILL  MOVE  mING  1/4  ChURD  TO  KING  HINGE  LINE 

"  =  i  PROGRAM  mILL  MOVE  RING  1//  CHORD  TO  mTn-J  HINGE  LINE 

•  x  h  INPUT  A*1NG  AS  STATION  UF  RING  TRAILING  EDGE. 

NTPOS  a  1  INPUT  Xl'AIL  AS  STATION  UF  TAIL  LEADING  EDGE 

■  =  2  PROGRAM  «1LL  MOVE  TAIL  1/4  CHORD  10  TAIL  tilNGt  LlNt 

•  s3  PROGRAM  MILL  MOVE  TAIL  1  /2  CHORD  TD  TAIL  b  1 NGE  LINE 

•  s  4  INPUT  XIA1L  AS  STATION  OF  TAIL  TRAILING  EDGE. 

HCDVT  *  1  NU  EFFECT. 

•  *  2  USE  NACA  KEPT.  1  2S3  ( HCDEV ITT )  FOR  FIN  ALJNfc  CLA. 

GOOD  FUR  RECTANGULAR  FINS  ONLY. 

GOOD  FUR  MACH  NUMBERS  BELOm  1.1  ONLY. 

NOUI  s  1  TRIM  OUTPUT  unLI. 

*  2  ALPHA  UUTPUT  UNLY. 
s  3  TRIM  AND  ALPHA  OU1PUI. 

NRUN  s  INITIAL  RUN  NUMBER.  -  I  TO  4  DIGITS. 

NDBASE  si  NO  EFFECT. 

-  2  OUTPUT  DATA  BASE  FILE. 

. .  INPUT  3 

NOSE  s  1  OGIVE  NUSE 

"  s  2  CONICAL  NOSE 

XLN  s  NUSE  LENGTH 

XL  s  TOTAL  LENGTH 

. . . . .  input  4 

RB*  s  RADIUS  OF  BODY  AT  KING 

RBT  s  RADIUS  OF  BODY  AT  TAIL 

DREF  s  REFERENCE  LENGTH 

Alter  •  REFERENCE  AREA 

DREFFT  s  REFERENCE  LENGTH  •  FT. 

AREFFT  s  REFERENCE  AREA  IN  SO.  FT. 

tit  *  TRANSVERSE  MOMENT  UF  INERTIA  -  SLUG-FT2. 

........ . . . ,,  INPUT  5 

NBOOY  ■  NUMBER  OF  POINTS  FUR  BODY  ALONE  CNACY  AND  XCPCY 

(IF  NBUOYSO,  OMIT  INPUTS  6.  PROGRAM  WILL  COMPUTE  SLOPE  DATA.  } 
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INPUTS  6 


I  •  A  C  h  *  TABLE  UF  MACH  NUMBERS  FOR  CNACY  AND  XCPCY 

IC*«ACi  =  1 A  Blit  OF  BODY  ALONE  NORMAL  FORCE  COEFF  SLOPES  IttASED  U« 
BjuY  CROSS  SECTIONAL  KtftKENCE  AREA)  -  1/KAD 
IXCPCJ  =  TABLt  OF  BODY  ALONt  CENTER  OF  PRESSURE  -  CAL  FROM  NOSE. 

. . . , . .  I  NPUl  7 

>.<sTu*0  PROGRAM  MILL  NOT  ADD  bOAfXAlL  CN A  AND  XCP. 

*  *1  PRUGk AN  USES  6UA1TAL  CARDS  TO  CALCULATE  BOAT  TAIL  ERECTS 

-BUD  =  DIAMETER  RAHu  (BOA  I T A 1  L/C Y L I NDtK ) 

XuBoL  a  RATIO  (BOAT  TAIL  LENGTH /CYLINDER  DIAMlTEk) 

. . . .  INPUT  A 

N*«F  s  NUMBER  OF  POINTS  EUR  XL'G  VS  MACH  nU«HFK 

. . . . .  inputs  9 

T*»F  =  TABLE  UF  MACH  NUMBERS  FOR  XCG 

1X0,  s  TABLt  UF  CENTER  UF  GRAVITY  -  CALIBERS  FROM  NOSE. 

I aElGH  *  TABLE  OF  *E1GHT  -  LBS. 


Input  id 


NFMA  s  number  OF  FREE  STREAM  MACH  NUMBERS 

. . .  INPUTS  11 

TFMA  *  TABLE  OF  FREE  STREAM  MACH  NUMfcERS 

TALI  *  ALTITUDE  FUR  EACH  MACH  NO,  -  t T  ABOVE  SEA  LEVEL 

THE  s  KEYNULDS  NUMBER  FUR  EACH  MACH  NO.  -  MILHUNS/FT. 

NOTES  PROGRAM  MULTIPLIES  IRE  BY  1  MILLION  TO  GET  RE. 

. . .  INPUTS  12-26 


TBU* 

s 

TABLE 

OF 

*ING 

ICR* 

a 

1  ABLE 

OF 

MING 

TTK* 

s 

TABLE 

01 

*ING 

TSaTE* 

a 

TABLE 

OF 

alNG 

TXNING 

s 

TABLE 

OF 

alNG 

THING* 

■ 

TABLE 

OF 

*InG 

TBOT 

a 

TABLE 

OF 

TAIL 

TCRT 

a 

TABLE 

OF 

TAIL 

TTRT 

M 

TABLE 

or 

TAIL 

TS*TET 

m 

TABLE 

OF 

TAIL 

TXTA1L 

a 

TABLE 

OF 

tail 

THlNGT 

a 

TABLE 

OF 

TAIL 

a 

a 

0.  (PROGRAM  SI 

TDELT 

a 

TABLE 

OF 

TAIL 

TDEL* 

a 

TABLE 

or 

*ING 

TALP 

a 

TABLE 

OF 

ANGLI 

EXPOSED  SEMI -SPANS (LOAD  0.  IF  NO  *InG) 

ROOT  CHORDS 

TAPER  R ATI OS ( CT/CK ) 

TRAILING  EDGE  SaEEP  ANGLES (DEG) 

STATIONS  (MEASURED  FROM  NOSE)  (SEE  NaPUS) 
HINGE  LINES 

EXPOSED  SEMI-SPANS (LOAD  0.  IF  NO  TAIL) 

ROUT  CHORDS 

TAPER  RATIOS  (CT/CR) 

TRAILING  EDGE  S*EEP  ANGLES  (DEG) 

POSITIONS  (MEASURED  FROM  NOSE) 

HINGE  LINES 


OF  ATTACK  -  DEG. 


. .  INPUTS  27*29  (USED  ONLY  IF  NUBASE*!) 

ICU)  a  FIRST  PACKED  WORD  FOR  COEFFICIENTS. 

IC(2)  a  SECOND  PACKED  aORD  FOR  COEFFICIENTS. 

IP  a  PACKED  WORD  FOR  PARAMETERS. 

(SEE  DATA  SAKE  INSTRUCTIONS  FOR  DEFINITION 
Of  PACKED  WORDS) 
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ul*fcN61UN  IMF  ( lb}  f  TXCGUS)  ,TmACH(15)  ,1CnACY  lib) ,  UCPCY  (lb) 
DiMENSIONTFMAUb)  ,  ICRM(lb)  ,  ilK"Uo)  ,TBU*(ib)  ,TS*rE**(lb)  ,lA*INGUb) 
DIMENSION  TCRf  (lb)  ,TTRT(lb)  ,1B0T(  lb)  ,TS*rEl(lb)  .TXTAlLUb) 

DIMENSION  THlNGR(lb) ,THlNGT(lb) 

DIMENSION  TREU6)  .TALTilb)  ,  TDELlUo)  ,TOEL«(  10)  ,  1  AbP(2u)  ,T«F.IGH(  1  b) 
DIMENSION  1F1LE(3) /TITLE  l  20) ,HUL110) 

DIMENSION  ACON(i) , R ( 32 ) , I ( 2b , 3 ) 

COMMON  /  bUDY  /  Xb.XbNjKE.DKEFFl ,OCK  , AC Y # A REF , NUSE , CN ab , ACPb 
COMMON  /MING/  bu*  ,  RBw  ,CK*  ,  X*J.NG,  XCEN*  ,  PHI*  ,  HInO  ,  A*AK  ,  AH* 

COMMON  /  4  i  N  G  2  /  CLA*,XCP*,bKAbfbNb*,5RW6,SKH* 

COMMON  /TAIL/  B0T,KBT,TR1 *XTAIL, XCENT.PHIT.HINGT, AT AH, ART 
COMMON  /TAIL  2/  CLAT,  XCP  T  ,  bR  T  t* ,  bK  h  f  ,SKiH,SRbT 
COMMON  /PARTS/  CNBY  ,XCPbl  ,  CN*  ,  XCP«  N  ,  CNT  ,  XCPTL  ,  Cl.  f  F 


REAL 

102 

103 

10b 

109 

110 
111 
112 

198 

199 

200 
202 
203 

209 

210 
211 
212 

213 

214 

215 

216 
217 
219 

219 

220 
221 
222 

223 

224 

225 

226 
227 

230 

231 

232 

2322 

2323 

233 

234 

235 

236 

301 

302 
400 


LAbEL 

FORMAT 

FORMAT 

FORMAT 

FORMAT 

FORMAT 

FORMAT 

FORMAT 

FORMAT 

FORMAT 

FORMAT 

FORMAT 

FORMAT 

FORMAI 

FORMAT 

FURMAT 

FORMAT 

FORMAT 

FORMAT 

FORMAT 

FORMAT 

FORMAT 

FORMAT 

FORMAT 

FORMAT 

FORMAT 

FORMAT 

FORMAT 

FORMAT 

FORMAT 

FORMAT 

FURMAT 

FORMAT 

FORMAT 

FORMAT 

FORMAT 

FORMAT 

FURMAT 

FORMAT 

FORMAT 

FORMAT 

FORMAT 


(2F10.4) 

(3F10.4) 

(BF10.4) 

(1312) 

(12. BX,3F10.4) 

120A4) 

(212, 6X, 3F1D.4) 

(1H  ) 

(1H1) 

( 1 H 1 , ' I NPU T  CABLES'/) 
(11  .2F10.4) 

(1M  ,  3 F 10*4) 
(1H0.20A4) 


(1H0, 
(1H  , 
(1H  , 

(1H  , 

(1H  , 

(1H  , 
(1H  , 
(1H  , 

(1H  , 
(1H  , 
(1H  , 
(1H  , 
UH  , 
(1H  , 
(1H  , 
(1H  , 
C1H  , 
(1H  , 
(1H  ,' 
(1H  ,  • 
UH 
(1H  , 
(1H  , 
UH 

UH  ,  • 
UH  ,' 
UH 


GEOMETER  1C  V  AH  I AbLE5 1 30X ' BOO Y ' bX ' *  1  ' bX 


CH 

CT 

bU 

AR 

TR 

S*bE 

S*Tt 

CBAR/CH 

YBAK/BO 

XM/CR 

CENT/CR 

RB 

XFIN 

XCENT 

A/AR 

HINGE 

PHI 


i*n  f  r n * 

UH  ,  ' XL 
(lH  ,'DREF 
UH  ,  'AREr 
(lH  , 'DREFFT 
(lH  ,'AREFFT 
UH  ,  •  Y1Y 
UH  ,'OBOD 
UH  ,'XLBOO 
(lH  ,'OCYODR 
UHO'PROBLEM 
( 1  HO  *  b ARMING t 


HOOT  CHORD* 

TIP  CHORD* 

EXPOSED  SEMI-SPAN' 

ASPECI  RATIO' 

TAPER  RATIO' 

LEADING  EDGE  S«tEP  ANGLE  (oeG) 
TRAILING  EDGE  5<*EEP  ANGLE  (CbG 
H  MEAN  GE0METER1C  CHORD/CR ' 

0  SPAN* ISE  LOCATION  OF  CBAK ' 
CHORD 4ISE  LOCATION  OF  CB AR ' 

R  DIST,  FROM  L.E.  CH  TO  CEnTHUID 
RADIUS  OF  BODY  AT  Fl *' 

DIS1.  FROM  NOSE  TO  L.E.  CH* 
DIST.  FROM  NOSE  TO  CENTRUIt) ' 
AREA  RATIO  (2  FINS/AREF)' 

DIST.  FROM  NOSE  TO  HINGE* 

FIN  ROLL  ANGLE  (DEG)' 

BODY  LENGTH* 

REFERENCE  LENGTH  • 

REFERENCE  AREA  • 

REFERENCE  LENGTH-FT.' 

REFERENCE  AREA  FT**2 ' 

MOMENT  OF  INERTIA' 

B0ATTA1L  D1A.  (CAL.)' 

BOATTAIL  LENGTH  (CAL.)' 

RATIO  (CYL.  Di A/REF •  OIA)' 

M1TH  INPUT  DATA  FILE-(TAPE5) ,  S T A 
TAPE5  NOT  CLOSED'/) 


, 'TAIL'/) 

3bX2F 1 u . 4 ) 
37X2F 10.4) 
2HX2F 1U.4J 
3  <*  A  2  F  l  0  .  ** ) 

3BX2F1U.4) 
•  lbX2F 1U.4) 
)'  lbX2F10.4) 

22A2F10.4) 
21X2 FID. 4) 
20X2F1U.  1) 
'  lbX2FlU.4) 
25X2F 10.4) 
20X2F1U.4) 
19X2F10.4) 
22X2 F 10,4) 
22X2F10.4) 
26X2F1U.4) 
24XF10.4) 
17XF1U.4) 
1UXF12.4) 
lbXF10.4) 
15AF10.4) 
lbXFlO.4) 
lbXF10.4) 
13XFJ0.4) 
10XF1O.4) 
rUS«' 14/) 


302  FORMAT  (  1  HO » MARKINGS  TAPES  NOT  CLOSED'/) 

400  FORMAT  (1H0, 1 1 X , ' BRUP ' , SX , ' bKCO ' , 5X , ' CLA ' , 5X , • CP • , BX , 

♦  'CNA' ,bX, 'XCP* ,5X, 'CMACG* , 5X, 'CMO* ,bX, 'CLP' ,bX, 'CLD' , 

♦  bX,'CMA',bX, 'CHD'.bX, *CN0',5X, 'CMDCG ' ) 

401  FORMAT  (1X#'B0DY  ',3bX,4F9.3) 
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402  FGMMAi  (lA,'wl*G  »,14F9.3) 

40  J  FORMAT  (IX,  '  T  A 1  Li  *  ,  1 4  F  9  «  3  ) 

404  y QRMAT  ( IX  # 'TOTAL* ,3bX,5F9.3) 

4444  FORMAT  (  1  X  ,  ' S IGA A  ',F9.3) 

404U  FORMAT  ("  •  ,T4J  ,  4»>  l  T4  3  T«8  , 

♦  T4j,"*-,T58, "FLIGHT  C0Uulfl0NS",TB8,"4«,/,l43,<,4",r88,"4" 

♦  ,/.T43,46("«<,n 

4041  FuRMAT  1 1  MO  ,  "  ",144, 4b  (  " 4"  )  ,  /  ,  1  4  J  ,  "  ♦  "  ,  T«  »  , 

♦  T4J,"*",T52,  "ANGLE  OF  ATI  ACK  DATA  (  SLOPES  )",  1  H  6  ,***,/,  143  ,  M  »  "  ,  lb  8  , 

♦  T4J,"4",  TB8,  *♦",/,  T43,4b("*" ) ) 

4042  FORMAT  l*  ",T43,4b("4"),/,T43,»«",T«8,**,\/,T4J,"*",T88,"4",/, 

♦  T4J,"»",Tol,"TKIM  PUlwl",  188,  T4J,  "4",  188,“*",/, 

♦  T43,**", T88, "♦",/, T4J,4b("*")) 

4043  FURMAT  11MO/  "  ,T44,4b("4-  J  ,/ ,T43,  "**,  [He  T43,  Irtb  ,  M*"  ,/, 

+  T4j,"*"f 155, "ANGLE  Of  ATTACK  S* EEP" , T«H f4 3 , f« a , 

♦  T43,"*",i88,"4",/,r43,4b{"*"); 

410  FURMAT  (1 HO, 'FMACH*' ,F9.i,2X, • ACC* * , F 9 . 3 , 2 A • * E luH {= ' , F9 . 3 , 

+  2  A  ,  •  AOT='  ,F9. J,2X,  'Kt*'  ,F 10.0) 

450  FURMAT  1 1H0,  JA,  •  ALPHA  »  ,4X,  'DLL*  '  ,4X  ,  *  t'fcLT  *  ,  bX  ,  •  C  4 •  ,  oX  , 

♦  1  ACP  1  ,  bX  ,  *  CmCG  1 ,5X,  '  HM*  •  ,oX,  'HMT'  ,5X,  'FTP  A  I  ’  ,  U,  'GLOAu*  ) 

451  FORMAT  11H  , 1 F8 . 2 , 3F9 . 3 , 1  OF  8 . 3 ) 

452  FORMAT  (iHO.'NU  TRIM  AFTER  1 , 1 3 , 2X , 'ITERATIONS  FOR  DfcL**', 

F7 , 3 , 3X  ,  ' DELI* '  ,  F  7 , 3  ) 

453  FORMAT  ( lho, 3X, 'ALPhA • ,bX, 'CN ', bX. 

+  'ACP' , bX , 'CMCG' ,bX, 'Ch* • ,5X, 'CHI' , 3X , 'FTRAT  '  ,4X, 

♦  'CNB'  ,4X,  '  ACPb'  ,5X,  'CN*  '  ,4X,  'XCP*'  ,5X,  'CM  '  ,4X,  'XCPT  '  ,  4A,  'Cl»Tll '  ) 

454  FORMAT  C1hO,'TK1M  PUlNl •  ) 

455  FORMAT  llH0,'NG  TRIM  HELD*  «LP=  ',  F7 . 3 , 2X  ,'  FOR  DEL,  *=  '  ,  F  /  .  3  , 

♦  2X, 'AND  DELT® • ,F7.3) 

45b  FURMAT  (1H0,'RUN  nUMBEk',14) 

990  FORMAl  (iHO.'UPE  INPUT  FILE  NAME  -  1  TO  1 2  CHARACTERS'; 

991  FORMAT  13A4} 

995  FURMAT  ( 1 H 1 , ' COnF IGOR A llON ' 1 4 ) 


,  1  »RITE  11,990) 

REAO  (1,991)  IFILE 
CALL  CLOSE  15, 1ST  A ) 

CALL  OPEN*  ( 5 , IF1 LE , 1,0,0,1S1A) 

IF  (ISTA.EQ.O)  GO  TO  8 
“RITE  (1,301)  1STA 
PAUSE 
GO  TO  1 
|  8  NCONF*0 
READ  (5,111)  TITLE 

READ  (5,4)  N*POS,NTPUS,MCDVT , NOUT , NRUN , NDBASE 
READ  (5,*)  NOSE , XLN , XL 

READ  (5,4)  K8m,RBT,DREF,AREF,DREFFT,AREFFT,¥1T 
READ  (5,4)  NBOOX 
IF  (NBQDY.EU.O)  GO  TO  4 

READ  (5,4)  (TM ACH (I),TCNaCY (i),TXCPCY (1 ) ,  lal,MBOOl) 
4  READ  (5,4)  NBTL,DBOD,XLBOD 
REAO  (5,*)  NMF 

REAO  (5,*)  (TMF(I),TXCG(I),T*E1GH(I),  1*1, NMF) 

REAO  (5,4)  NFMA 

REAO  (5,4)  (TFMA(l),TALT(l),TRb(I) ,  1*1, NFMA) 

CALL  READX  ( 5 , 80 , HOL , TBO* , NH , NBO* , EOFF ) 

CALL  READX  (5,80,HOL,TCR«,Nh,hCR*,EOFF) 

CALL  READX  (5,90,HOL,TTKM,nH,NTRn,EOFF) 

CALL  READX  (5 , 80, HOL, TS*TE» , NH , N5MW,E0FF ) 

CALL  READX  (5,80,H0L,TX*lNG,NH,NX*ING,fcUFF) 

CALL  READX  (5,80, HOL, THING*, NH,NHINGM, EOFF) 

CALL  READX  ( 5, 80 , HOL , TBOT, NH , NBOT , EOFF ) 

CALL  READX  (5, 80,H0L,TCRT , NH , NCRT , EOFF ) 
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CALL  kL  ADA  (5,6o,HOL,TlHl  ^H.NTKT.EOf'F) 

CALL  RtADX  (5,80,HUL,rSwrtT,*B,NSwT,LUH' > 

CALL  READX  (6,«0,HUL,  fXT A 1 L , , N Al A 1 L , EUFF ) 

CALL  KKADX  (5,80,H0L,THING1  ,  NH , NH1 NGT , LUFF ) 

CALL  HEAOX  (5,80,HUL,T0fcLT,Nb,NDtLT,tUH- ) 

CALL  hfcADX  (5,60,HUL,TDfcLw,NM,NDELw,LUH  ) 

CALL  rtADX  (6,bO,HOL,  TALP  ,  NH  , ALP  ,  EUF t  ) 

IF  (NjbASE.tU.2i  CALL  DA1A1 
CALL  CLOSE  ( 5  » 1ST  A ) 

IF  (lSTA.Nt.D)  "MU  (1.302) 

. .  wKITt  INFO  I  T  AbLES 

WHITE  (6,200) 

IFCNbULl .EU.O)  GU  TO  2 

WHITE  (6,203)  (TWACH(l),fCNACK(l)#rxCFCV(I),  ial.NBuUl) 

2  wHITt  (6,198) 

WHITE  (6,203)  ( IMF U ) ,TXCG( I ) , TwtlGH ( I ) ,  1*1, NMM 


PI«J. 14159 
RAD*57. 295)8 
SMALL* 1 . E“ 30 
PHlNSO.O 
PmITsO.O 
DCI*2.»RBT 

dctodr*dcy/dref 

ACX*(P1/4)4DC**»2 


i 


i 


00  90  12*1, NCRW 

CKw*TCRw (12) 

DO  9o  13*1, NTRW 

THw*TTRw (13) 

DO  90  14*1, NBOW 

BUW*Tt)Uw  ( 14  ) 

DO  90  15*1 , NSWW 

SwTEw*TSwTEw(15) 

DU  90  16*1  ,  NXwlNG 

XwlNG*TXwING(I6) 

DO  90  112*1, NHINGw 

HINGw*THINGn(I12) 

* 


CALCULATE  WING  GEOMbTKl 


CALL  ROLL  (BOw,CRW,TKW,RBW,DKEF,ARLF,CLPWl ,CLD*1 ) 

CALL  FINGEOM  (CRw , TRW , 80w , AREF , SwTEw , CTw , AHW , SwLEW , CB AHw , T b AH « 
♦  XMw, CENT*, AFAR, TANLw.SwHAw) 

AwAR«2  ♦AFAR 

IF ( NWPOS , EO . 2 )  XwlNG*HINGw-( .25*CBARwtXMW)*CRN 
1F(NwP0S.E0. 3)  XWlNG*HINGw-( ,504CBARR4XMw)4CRw 
IF (NwpUS.EO.4)  XwING*XwING  -  CRM 
XCCNws(CENTw*CR*+X* 1NG) 


DO  90  17*1, NCRT 
CRT*TCRT (17) 

DO  90  18*1, NTRT 
TRT»TTHT(I8) 

00  90  19*1 »NBOT 
S0T*TB0T(I9) 

DO  90  110*1, NSMT 
SMfET*TSMTET(IlO) 

DO  90  111*1, NXTAIL 
XTAIL*TXTA1L( 1 1 1 ) 
00  90  113*1, AHINGT 
HIMGT*THINGT  (Hi) 


CALCULATE  TAIL  GEOMETRT 
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CALL  ROLL  (BUT, CRT, Ikl ,RB'i ,OKfcF , AHfcK ,CLPTi ,CLuT 1) 

CALL  FihGEUM  (CRT,  f  K 1  ,  bU  l ,  AKEF  ,  S*TtI  ,  Cl  J  ,  AH  f  ,  S»»LE.  i  ,  Ctt  AH  1  ,  t  B  AK T  , 
♦  XMT.CtMT,  AF AR  ,TAnl1,SwHAT) 

ATAR*2.*AFAk 

IF ( MTPOS , EU . 2  )  XIAIL*HINGI-( . 25*C& ART* AMT ) *CKT 
IF(NTPOS.tu.i)  ATAlL*HlNGl-(  .604CBART  +  XM1  )*C«1’ 

IF(hTPuS.tU.4)  ATA1LSA1A1L  -  CRT 
XCENT*ICEMTWCRI*aIAIL) 

lF(HXNUr.bU.O.>  HINGTaXlAiL  *  (AMI  ♦  .  2  5  *  C  k  A  k  T )  *  C  R I 

. . .  GtGftik*  output 

NCONF  SNCON*  *  l 
WHITE  (6,995)  WCONF 
WHITE  («,2U9)  TITLE 
*RI  It  (o,210) 

-RITE  (6,230)  Xu 
wRITt  le,23l)  Ukt* 

WRITE  (6,232)  Akt-F 
wRITF  (6,2322)  ORt.FFT 
WRITE  (6,2323)  AREFFT 
WRITE  (6,233)  Hi 
IF  (N6TL.E0.0)  GO  TO  5 
WRITE  (0,234)  OdOD 
WRITt  (6,235)  XLBOU 
WRITE  (6,236)  DCXODK 

5  -RITE  (6,213)  1 BU« ( l 4 ) , TBO l( 19 ) 

WRITE  (6,211)  TCh»(I2),TCRT(I7) 

WHITE  (6,212)  CTw,CIT 
WHITE  (6,215)  Tkw  , ThT 
WRITE  (6,214)  A  R  » , A  R 1 
WRITE  (o,216)  Sw  LE* , S w  L 1 1 
WRITE  (6,217)  SwTEw.Swrtl 
WRITE  ( 6 , 2 1 B )  CBAkw , CB ART 
WRITE  (6,219)  1 b Ah  w , j  BART 
WHITE  (6,220)  AM w,AMT 
WRITE  (6,221)  CENT*, CENT! 

1,  WRITE  (6,222)  Raw, RBI 
i  WRITE  (6,223)  Awing, X1A1L 
WHITE  (6,224)  XCfcNw,xCENT 
WHITE  (6,225)  AwAR, ATAR 
WRITE  (6,226)  MlNGw, HlNGi 
WHITE  (6,227)  PHI w , PHIT 

..................... . .  START  MACH  VUMBER  uUOP 

DO  91  11*1, NFRA 
FMACH*TF»*A(I1) 

ALT»TALT(11 ) 

RE*TRE ( 1 1 ) *  1 , Eb 

CALL  STAN62  ( ALT ,CS, RHU,PHESS, TEMP , 2 ) 

0«( .7*PRESS*FmACH**2) 

BETA*SURT (ABS(FMACH**2*1. ) ) 

CALL  IhTERP  (2 . 1 , XCG , TXCG , FM ACH , IMF , NMF , 2 , m IN , * A A ) 

CALL  IHTERP  (2»l, WEIGHT, TwE1GH,FMACH,TMF,NMF, 2, mIm, MAX) 
XCG*ACG*DCX 

IF  (WEIGHT. EG. 0. )  WEIGHT*!. 


START  ZERO  ALPHA  LOUP  •  SLOPES 

. . . . . . .  HOOT  FORCE  COEFF. 

IF  (NBODY.NE.O)  GO  TO  7 
XLA*XL  -  XLN 
XLN2*XLN/DCI 
XLA2«XLA/0CK 

CALL  BOOT  ( NOSE, X LN 2, XL A 2, F MACH, CM  ACT  ,  XCPCY ) 
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GO  1U  3 

7  CALL  WfckP  W,  1  ,C"«ACK  ,  iCNACK  ,*  MACH,  lMACH  ,  UtvJul  ,  S  ,  -11  n  ,'4AA  J 
CALL  IiiTk.-r  (2, 1 , XCPCI,  l  ACPC*  ,  t  M  ACH  ,  TNACH  ,  NbUOK  ,  3,  Hl.»  ,  4  AX  j 
3  cu'-ruot 

XCPC*  =  XCK l *OC< 

CNACKsCNACl*ACi/ ARfcK 
CNABTsO. 

XCP  BT*<». 

IF  (NttTL.iu.O)  GO  i'U  6 

CALL  ts 0 A T  L  (MACh  ,  DtSOU,  ALoLO,CNAaT  ,  XCPttT  ) 

XCPB'isACPb  r*l>C<UOH*UHLF  ♦  AL  -  XLBOD*PC  m>K*L>REF 
CNABr*Cr»Abl*AC</AK£F 
6  CNAe*CNACl  ♦  C^ABi 

CMA8s(CtfAC**ACFC](  4  C*Ad 1  * ACPO 1 J / LREF 
XCPB=C*Ab«t)KLF /(CNAB  4  S*AbL> 

AhHCXs ( XCG<*XCPC i ) /  UREF 
ARMBls(XCG-XCPBl)/uR£F 

CMUBs-2.*(CNACf *ak*C***2  4  CNArtlAAHKtHl**^) 

CwACGB»CSAC**AKaC I  *  CNAbT*AFi*BT 

. . . .  *I!G  FORCE  COLFF  . 

F  w*Q  . 

XCP«*=U. 

CMUweO. 

CbaCGasu, 

CLPnsO . 

1F(BU*.EO.O.O)  GO  ro  561 

CALL  BKIT  (FMACH,  AR« , TR* , S*HA » , CLA N , CPC» ) 

CPw*CPCB«CbAR«  ♦  AA* 

IFUACDVT.tU.2J  .AND.  (FMACH, LE.  I  .  1  }  )  CALL  MCOEVT  l  F*  ACrt  ,  AH  *  ,  CLA  a  ) 
CLPa*CLP«J*CLA«/(PI*AKi*/2.  ) 

CLD**CLD*l*CLA*/(PI*AR«/2. j 
XCP*«(CFa*CH»  ♦  XAJNGJ 
BM*s99.9999fc4>3U 

IF  (TANLU.GT.O.D)  B**a0tTA/l'A*L» 

CALL  1NTSLB  (  1  .MACH.BtlA,  8R*  ,  k  6  A  ,  Ck  A  ,  1  rt  *  ,  bO»  ,  AH  * ,  Cl  a  *  , 

1  BKNBfbAB*.SA*BtSKB*J 
F  **  (  BKNB  ♦  BKB*)*CLA*»A*iAH 
ARNNH*(HlNG»-XCp..)/DREF 
CHA«fsBKwto*CLA«f*(  AaAR/2.  J  *ARN"b 
C:!D*=i>K«B*CLA*C*iAi*AR/2,  )*A R*«H 
CNOns(SMB*SKB»  J*CLAk*A«AR 
AKMN*( XCG"XCP* ) /DREF 

cmdhcg*Cndn*ara« 

CMOnk*2.«(F«  *ARN»i**2) 

CMACG»*FB  *ARM*> 

5b  1  CONTINUE 

. . . . .  TAIL  FORCE  COEFF. 

FTI«0. 

XCPT*U . 

CNOT*0, 

CMACGTsO, 

CLPT*0. 

irCBQT.EU.O.O)  GO  TO  651 

CALL  BRIT  (FNACH, ART .TRT.SNHAT 'CLAT.CPCb ) 

CPTaCPCB*CBART  ♦  XAT 

IF ( (MCDVT , EU, 2) . AND, l FMACH . Lt .1.1))  CALL  KCDEVT  (FA ACH. ART, CL AT) 
CLPT»CLPT1 *CLAI/ ( PI * ART /2 . ) 

CLDT*CLOTl*CLAl/( PI* ART/2. ) 

XCPX«(CPT*CRT  ♦  KTAIL) 

8RT«99.9999E430 

IFCTANLT.GT.O.OJ  B*T«BETA/TANLT 

CALL  1NTSLB  ( 1 ,FNACH, BETA ,  BAT , RBT ,CRT ,TRT ,BOT , ARI ,CLAT , 

I  BKTB'BKBT .SKTB'SKBT) 
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dIAlLsXL-AlAIL-CKT 

XXLI»*ALbGD*ORtfr 

DDBaLBOMUKfc) 

CAUL  AhAKF  MACH,  CRT,  BrAlL,XXLB,DL)b,  OCT  ,A1A) 

BKB1*BKbT  *A 1  a 
5KBTsSAo1*AIA 

FTa(bK  lb  ♦  bKBn*CLAT*ArA« 

CNl)T*(SKIb+SN8T  )*CLAT*A  T«K 
CADTCGaCNOI*(XCG-ACPl )/UKtF 

. . . . . .  xuRTtX  CALCULATION 

CLI Asu , L 
Xlao.O 
FwS»*0 , 0 

ARMrs(  XCG-xCPI )/DkEF 
ARMTHaCHlNGi-ACP  r )  /DKEF 
IF(bGw.EU.O.O)  GO  TO  655 
XHwatXCPwXwlNG 
XBART*-»XCPI-XfAlU 
ALP*. 5 

CALL  FACTOR  (bOw  ,  Rbw ,  CR  w  ,  XH  A  ,  X  w  1  N(i  ,  o  .  ,buT,KBT  ,  T  K  T  , XT  H I L ,  XtJART  , 
1  ALP, XI, MS*) 

FwMRw  s  MS**BU» 

CL1A*ICLA**CLA1 *BMB*Xl WhUT /(2.*PIAARI*R WMRw))*AwAK 

CLlOASCCLAd^CLATASKbBOXlAbOl'/l^.^PlAARTAFwMK*)  )*X*AR 

CMIDwC*CLiU**ARMT 

b55  CONTINUE 

FTI*H  ♦  CLI  A 

FTKAT*f TI/IFI+SMALL) 

CHAT*(bM'B*CLAX*ATAR  ♦  CLlA)*.5*AHwTH 
CnOTe(SKTb*CLAT*ATAR  ♦  Cl  1 A  )  *  . b*  ARmTh 
1  CMt)Ta-2.*(FTI*AKMT**2  ) 

:  CMACGIaFTl*ARMX 

651  CONTINUE 

|  IFUbOT.EU.O.)  .UR.  IBOm.EO.O.  )  )  GU  TU  b52 
I  CND**CNDw*CLID* 
j  CMUwCGaCMDwCG+CMlDwC 

652  CONTINUE 


IF  (MJUT.NE.t)  WRITE  (6,199) 

WRITE  (b,  198) 

WRITE  (6,199) 

WRITE  (6,4040) 

WRITE  (6,410)  FMACH,XCG, WEIGHT, ALT, RE 
WRITE  (6,4041) 

WRITE  (6,400) 

WRITE  (6,401)  CNAB,XCPB,CMACGb,C*OB 
IF  (BOw.EQ.O)  GO  TO  653 

WRITE  (6,402)  BAwB,BKB»,CLAw,CPw,Fw,XCPw,CWACGw,CAQw,CLPw, 
♦CLOW, CHAW, CHOw, CNOw.CMOwCG 

653  IF  (BOT.EQ.O)  GO  10  654 

WRITE  (6,403)  BKTB, BK8T , CL AT , CPT , F TI , XCPT , CMACGT , CMUf , CLPT , 
♦CLOT, CHAT ,CHOT ,CNDT, CM OTCG 

654  CONTINUE 


COMPLETE  CONFIGURATION 
CNAaCN AB  ♦  Fw  ♦  FTl 

CNA*CMAB  ♦  Fw*XCPw/OREF  ♦  FTI *XCPT/DRE) 
XCP«CMA*OREF/(CNA  ♦  SMALL) 

CMQbCMOB  ♦  CMOW  ♦  CMUT 

CMACG*CM ACGB  ♦  CMACGW  ♦  CMACGT  82 


CLP*CLPm  ♦  CLPf 

SlGMA*SuKl  (  AeSlCMACG)4u*AREFF  T*UK£F F1V  (I  I* ♦SMALL)  >  /  t  2  . ■* t> i  ) 

•Hilt  (6/4041  CNA,XCP,CMACG,O0,CLP 
•HITE  (6/4444)  SIGMA 


END  2EHU  ALPrtA  LOUP 


Jb  92  Ll*l / NUELT 
■jELTsTLELT  (LI  ) 

Ij  92  L2*l  /  NOEL* 

-tL^sTDEL*  ( L2 ) 

. . .  TRIM  CUNUH  l  JhtS 

IF  (NOUT.EU.2)  GO  TO  94 
4iiP*-(CMD*CG*DtL*  ♦  CM0TCG*0tuT )/CMACG 
STRIMalOO 
00  95  1TR*  1 / NTKlM 

CALL  HIGHALP  (FMACH,XCG,AlP, DEL*, PELT, CN/XCP/CMCG/CHa.CWi, 

♦  FTKAT,DKEF) 

UALPS*0 . 5 

IF  (ALP.LT.O)  DALPS»-0.5 
ALPS*ALP  -  DALPS 

CALL  HIGHALP  ( FMACH , XCG , ALPS , DEL* , CELT , XI , X2 , OCG l . X 3 , X4 . X5 , UKEF ) 
CNCGS* (CMCC*CMCG1 )/0ALPS 
OALP»CMCG/(CMCGS  4  SMALL) 

IF  (ABS(DALP).LE.O.l J  GO  TU  96 
ALPsALP  -  OALP 

95  CONTINUE 

•RITE  (6/452)  NTR I M , DEL • , UEL J 
GO  TO  94 

96  ALPMAXS4G. 

IF  (ABS(ALP).LT.ALPMAX)  GO  TU  97 
•HITE  (6/455)  ALPMAX , DEL* / DELT 
GO  TO  94 

97  GLOAQ*(CN*COS(ALP/RAu) )*9*AREFFT/*E1GHT 
rtM*»CHM*Q*ARtPFT*DR£FFT 
nMT»CHT*o*AHEFFT*DREFfT 

•RITE  (6/196) 

•RITE  (6/4042) 

•RITE  (6,454) 

•HITE  (6,450) 

•RITE  (6,451)  ALP, DELM/OELT/CN/XCP/CMCG/HMW/HMT,FTRAI/G LOAD 

. . .  ALPHA  CONDITIONS 

94  IF  (HOUT.EO.l)  GO  TO  92 

ACUN ( I  )  * '  • 

ACON ( 2 )* 1  * 

ACON(3)»»  • 

H(1)*NALP 
H(2)*NRUN 
H(3)«l 
H(4)»PHIT 
H($)»FRACH 
M(6)«RE/l,E6 
l«IK  (6,4043) 

M«in  (6,436)  NRUN 
MftITC  (6,433) 
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DU  9.)  i-  3  =  1  *  N  A LP 
ALP*1 AuPlL3  J 

CALL  hIGHALP  ( FMACH , XCG . ALP , DLL" , DEL I , Cw , XCP , CMCG , CH * , CRT , 

♦  FlrtAT,QH£t- } 

-RITE  t  b , 45 1  )  ALP,CN*XCP»CmCG,CHw*CHT*P  i  R  AT  , 

♦  Ob! ,XCPB* ,CN»,XCP*N,CNl ,XCP1L,CLTV 
UL3,  l)sALP 

UL3*2)xCn 
UL3,3)sCMCG 
93  CONTINUE 

IF  (NOBASE  •  £0.  2  )  CALL  DATA  (AC0N*H*Y) 

nRUN  *  NRUN  ♦  1 

92  CONTINUE 

91  Continue 

90  CONTINUE 

STOP 

END 

SUBROUTINE  HIGHALP  l  FMACh  ,  XCG  ,  ALPHA  ,  DEL** ,  L>E0  T  ,  CN  »  *CP  ,  C  4CG  , 

♦  CH«,CHT*FTRAI,OREFj 
HIGH  ANGLE  OF  ATTACK  AERO 

ALL  GEOMETRIC  VARIABLES  MUST  HAVE  SAME  SYSTEM  OF  JNlTS. 

I  ,£«  *  In,*  FT,*  CAL,*  ETC, 

COMMON  /BODY/  XL,XN,RE,0,DCY,ACY*AREF,N0SE,CNAH,XCPB 
common  /nING/  B0n,RB*,CR*i,X«1NG*XCENn,PH1N,HING*»,A*AR,AR*i 
COMMON  /«ING2/  CLAM*XCM*BK«b*6Kbw*SKN»,SKB« 

COMMON  /TAIL/  bOT,RBT,TKT,XlAIL*XCENT,PHIT,HlNGT, ATAR, ART 
COMhuN  / TAIL2/  CLAT,XCT,BKTb*8KBT,SKTb,SKBT 
COMMON  /PARTS/  CNHY,XCPB!,CNn,XCPm,CM,XCPT,CLI 

. . .  BODY  FORCES 


SMALL*!. 0E-30 
RAO*57. 2957b 

ALP*AbS(ALPHAJ 

XL2SXL/DCY 

Xn2»XN/DCY 

,XCG2*XCG/DCY 

CALL  CROSFLU  l Fm ACH , ALP , RE , D , XL2 , XN 2 , NOSE , XCG2 , CN V , CM V , X3 ) 
CALL  CROSFLI  IFMACH *  ALP *  1 , 1 , an2 , XL2 , XCG2 , X 1 , X2 , XCP V 2 ) 
XCPVSXCPV24DCY 
CnV*CNV4ACY/AREF 

CNAbD*CN AH/RAD 
CNB**CNABD4ALP  ♦  CN  V 
XCPBYSXCPB 

IFtALPHA, EG , 0 , )  GO  TO  44 

XCPBY*(XCPB*CNABD*ALP4XCPV*CNV )/(CNBYmSMALL) 

IFIALPHA.lt, U. )  CNBYs-CNbY 

. .  MING  FORCES 

44  IF ( BOM , NE,0 , 0 )  GO  TO  750 

CLMCO.O 
XCPtsO ,0 
CHMsu.O 
GO  TO  751 
750  CONTINUE 
Di«ALPHA  4  DEL* 

CLAMOaCLAM/RAD 
CL*A«CLANO*ALPHA 
CL«D«CbAHD*D£bR 
XCMaXCN 

ARMHaa ( HI NGN-XCPM ) /OREF 
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CL»*((  BK*B4  BKBN )  *Cl>« A  ♦  l  SKnB4  SKbN)*CLbD)*AwAK 
CH*a(BK«b4CL*A  4  SK*b*CL«*0)  ♦AkAH«i*A»A«/2  . 

CHM*CHh* ( CUS ID1/RAD)  4  S1N(D1/KAD)4*2) 

. . . .  TAIL  FORCES 

751  IF ( BUT .ftfc.O.O)  GU  XU  bbu 

CL  1 1=1) .  0 

XCPI*0.0 

CHT»0. 0 

GO  TO  651 

650  CONTINUE 
02*ALPh A  4  DELI 
CLATOaCLAT/RAD 
CLTA=CLATD»ALPHA 
CLTD=CLATO*OELI 
XCPtsXCT 

CLT*((  bK Xb4  bKBX)*CLTA  4  (  SKTB*  SKbT ) *CL1 U ) * A1 AH 

AHi**MlNG««Al»ING 

XbANT=4XCPT"XTAIL 

CLI*0, 

IF (bOw .Eu.Q. )  GU  TO  60 

CALL  VORTEX  (BO* , RB* , CR* , AH* , XN1NG , DEL* *  PH I  * , A* AK , BOX , RBI , IRT  , 

C  XTAIL,ABART, PHIT,  ART,  ALPHA, CLA* ,CLAT,BK*B,SK*B,CLI ,F*S* ) 

68  CONTINUE 
ARMHT*(H1NGT-XCPT)/DREF 

CHI=(BKT8*CLTA  4  SKTB*CLTO  4  CLI)*ARMHT*ATAR/2. 

CHT*CHT*(CUS(D2/RAD)  4  SIN(02/RAD)»*2) 

CLTI*CLT4CLI 

FTRAT«CLTI/ICLI  4  SMALL i 

651  CONTINUE 

. . . . . . .  COMPLETE  CDNFIGUHAXIU,. 

CnnsCLw 

C*T*CLT 

C*tl*CLTI 

j  CN  SCNBY4CNH4CNTI 

CM  =4CNBY*XCPBY  ♦CN*4XCP*4CNTIAACPT 
'  XCP*CM/(CN  4  SMALL) 

CMCG*(4CNBY*(XCG-XCPBY)4CN**(XCG-XCP*)4CNIH‘(XCG-ACPl) )/DREF 
RETURN 

Enu 

’  SUBROUTINE  DATA1 

!  INITIALIZE  DATA  BASE  FILE  INPUT. 

'N2  a  NO.  OF  COEFFICIENTS, 

N3  a  NO.  OF  PARAMETERS. 

CHARACTER*^  OFILE 

COMMON  /PACK/  XC( 2) , IP, N2 , N3 

200  FORMAT  (1H0, 'ENTER:  IC(1)  -  8  HEX.  CHAR,,  I.E.,  ABCUEF12 ' ) 

201  FORMAT  C1H0, 'ENTER:  IC12)  -  8  HEX.  CHAR.,  I.E.,  C0280000') 

202  FORMAT  ClHO, 'ENTER:  IP  •  8  HEX.  CHAR,,  I.E.,  AF2E1000') 

210  FORMAT  (28) 

300  FORMATC1HO, 'ENTER  DATA  BASE  FILE  NAME  -  I.E.,  DATA: 2CQPHED.DAT' ) 

400  FORMAT(AIB) 

401  FORMAT  ( lHO , ' PROBLEM  WITH  FI L£a * , A I  8 , 3X , » IS TATa • , 14 ) 


13  NRITEd  *  300) 

RE AD (1,400) OFILE 

OPEN ( 10, FI LEaQFlLC, FORM* 'BINARY ' , COUNTS )■ ' RECORD ' 

*, RECL*32, BLOCKS! ZE*32, I OSTATsl ST AT, ACCESS* ' DIRECT • 
4, STATUS* ' NEN • ) 

IP  (18TAT  *EQ.  0)60  TO  12 
NRZTE  (1,401)  OFXLC, X3TAT 
PAUSE  'CORRECT  PROtLEM  AND  CONTINUE* 
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Go  IU  13 


12  CONTINUE 

N2*J 

N3*6 


•HITE  (l,20OJ 
READ  lb, 210)  1C11) 

•RITE  (1,201) 

READ  (5,210)  1C ( 2 ) 

WRITE  (1,202) 

READ  (5,210)  IP 

RETURN 

END 

SUbSUUT I NE  DATA  (ACON,h,K) 

DIMENSION  1(25,3),  ACON(i),  H(32),  C(2,t>4) 

COMMON  /p ACM/  1 C(2),IP,N2,N3 
H  *  VALUE  OF  HEADING  PARAMETERS. 

Y  *  CGEFF.  ARRAY  AT  EACH  DATA  POInT. 

ACON  s  CONFIGURATION  NAME  -  3  WORDS  -  12  CHARACTERS. 

C  a  MAX  AND  MIN  VALUES,  (LOAD  ZERO) 

CORFILE.F1N  PROGRAM  DETERMINES  MAX'S  AND  MIN'S. 

DO  10  lxl,2 
DO  10  J*l,b4 
10  C(I,J)*0. 

WRITE  (10)  N2,NJ,1C,1P, ACON 
WRITE  (10)  ( H ( 1 ) ,  1*1, N3) 

WRITE  (10)  (C ( 1 , I ) , 1*1 , N2 I 
WRITE  (10)  (C(2,I),I*1,N2) 

NPT5*« ( 1 ) 

!  DO  50  1*1, NPTS 

150  *RITE  (10)  (T(I,J),J*1,N2) 

!  RETURN 
I  END 

SUBROUTINE  BODY ( NOSE , XLN , XLA ,* MACH , CNA , XCP ) 

INPUTS 

NOSE  *  1  OGIVE  NOSE 

NOSE  *  2  CONICAL  NOSE 

XLN  *  NOSE  LENGTH  (CAL) 

XLA  *  AFTERBODY  LENGTH  (CAL) 

OUTPUT 

CNA  *  NORMAL  FORCE  CUEFF,  (1/RAD) 

XCP  *  CENTER  OF  PRESSURE  (CAL  FROM  NOSE) 

TABLE  1  *  BODY  DATA  FOR  TANGENT  OGIVE  NOSE  WITH  VARIOUS  NOSE  LENGTnS 

TABLE  2  *  BODY  DATA  FOR  TANGENT  OGIVE  NUSE  WITH  VARIOUS  AFTERBODY  LENGTHS 

TABLE  3  *  BODY  DATA  FOR  CONICAL  NOSES  OF  VARIOUS  ANGLES 

TABLE  4  *  BODY  DATA  FOR  CONICAL  NOSE  WITH  VARIOUS  AFTERBOUY  LENGTHS 

DIMENSION  TMACH(15),TCNA1(15,3),TXCP1U5,3),TXLN(3),TCNA2(15,3) 

$ ,TXCP2( 15 , 3) ,TXLA(3),TANG(4) , TALA (5) , TCNA3(15,4) , TXCP3( 15, 4) , 

8TCNA4( 15,5),TXCP4(15,5) 

DATA  UMACH(l), 1*1, 15)/. 8,. 9,1., 1.05, 1.1, 1.2, 1.4, 1.6, 1.9, 2., 2. 5, 

•3. , 3.5,4, ,4.5/ 

DATA  (TXLN(X) ,1*1 ,3)/3.,4.,5./ 

DATA(TXLAd)  ,I*l,3)/4.,6.,9,/ 

OATA(TANGd)  ,1*1, 4) /6.  ,7. 125,8.  ,10./ 

OATA(TALA(l ) ,1*1,5)/4.,6.,8.,9.5,11./ 

DATA  ( (TCNAf (I , J) , 1*1 , 15) ,J*1 , 3J/.042 , .043, . 0465, . 048 , .0465 , .0455 , 

8. 0435, .0432, .044, .045, .05,. 054,. 055,, 0575,. 0575,. 043,. 044, .0465, 

8. 048,. 0445,, 044,. 042, *0412,. 0415,. 0425,. 047, ,0518,, 055,. 055,. 0545, 

8.044, .045, .046,. 048,, 0465,. 043,. 042,, 0405,, 0408, .041,, 0445,. 05,. 05 
83, .053, .054/ 
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DATAUTACPl(l,J),l*l,lS),0*l,3)/.65,.55,.25,.l7,.25,.b5,.97,1.06, 
$1*08, 1,05, ,66,, 7, ,5, ,3, ,06, 1,17, 1,05, ,55, ,45, ,65,1,, 1,4, 1,57, l9b, 
*1.5*, 1.46, 1.3, 1.1,. 9,  .65  ,1.6, 1,35, .9, .77, .96, 1,51, 1.93, 2. 1,2. 12, 
*2. 1,2. ,1.9, 1.69, 1.47, 1.2/ 

DATA((TCnA2(1,J) , I*1,15),J*1 ,37/0. 041 5,0. 044,0.0465,0. 048, 0.0465, 
$0.0445, 0.042, 0.041 8,0. 042, 0 . U425 , 0 . 047u , . 051 5, .054, .055, .056, .043, 
$0.045, 0.04b, 0.0505, 0.0495, 0,04 7, 0.044 6,0.044,0. 0445, 0.045, 0.049, 
$0.053, 0.0545, 0.056, 0.0575, 0.045, 0.047, 0.051, 0.052, 0.0515, 0.049b, 
$0.047, 0.0462, 0.046, 0.0465, 0.0495, 0.05 3, 0.0552,0. 05 7, 0.056/ 
DAtAllTXCP2U,J),  1*1,151,0*1,3 1/1.62, 1.4, 1.05,0.92, 1.1 5, 1.46,1. 95, 
$2. 1,2. 1,2. 08, 1.97, 1,8 1,1. 57, 1.3, 1.2, 1.1 5, 0.90, 0.40, 0.25, 0.045,1. U5 
4,1.56, 1.84, 1.90, 1.90, 1.76, 1.65, 1.4, 1.13, 0.96, 0.7^, 0.41, -0.10, -0.25 
$,0.0,0.75,1.4,1.65,1.70,1.70,1.55,1.4,1.2,1.0,0.62/ 

DAXAUTOA3U, J)  ,  1*1,15)  ,  J*1 , 4) /0.0405 , 0 .042 ,0 .046 , 0 . 04  75 ,0.0455, 
$0.04  4, 0.042, 0.04  1,0. 04 1,0. 042, 0.0  44, 0.046, 0.05  2, 0.0  53  5, 0.05  3, 0.0  40 
$5, .042, 0.046, 0.04 7 5, 0.0455, 0.044, 0.0 42, 0.0 4 2, 0.0 42, 0.0 4 3, 0.046, 
$0.05, 0.05 3, 0.055, 0.056,0. 040 5, 0.042,0. 046, 0.04 7 5, 0,0 455,0.044,0.0 4 
$4, 0.042, 0.043, 0.044, 0.047, 0.051, 0,0545, 0.057, 0.056, 0.0405, 0.042, 
$0.046,0.0475,0.0455,0.044,0.0436,0.0435 ,0.044,0.045,0.049,0.053, 
$0.057,0.059,0.060/ 

0A1AUTXCP3U,  J)  ,1*1,15)  ,  J*1 ,4) /0. 7  5,0,552, 0.20, 0.05, 0.325,0. 65, 
$0.94,1.05,1.1,1.09,0.90,0.65,0.42,0.24,0.20,0.55,0.40,0.05,-0.10, 
$0.15,0. 47, 0.70, 0.85, 0.85, 0.6 3,0. 60,0. 37, 0.15, -0.05, -0.11 ,0.37,0.20 
$,-0.10, -0.20,0. 0,0.30, 0.57,0. 7 0,0.70,0. 65,0.45, 0.22,-0.01 ,-0.20, 
$-0.27, 0.20, 0.05, -0,30, -0.36, -0.25, 0.05, 0.30,0. 45,0.4  7 ,0.4  3,0.242, 
$-0.05,-0.30,-0.48,-0.50/ 

0A1A (( TCNA4 ( 1, J), 1*1,1 5 ),J*1, 5) /0. 039, 0.041, 0.0445, 0.045, 0.0435, 

$0. 043, 0. 041, 0. 040, 0. 041 5, 0. 042, 0. 04b, 0.049, 0.051 5, u. 05 3, 0.054, 
$0.0405,0.042,0.046,0.0475,0.0455,0.044,0.042,0.042,0.042,0.043, 
$0.046, 0.050, 0.053, 0.055, 0.056, 0.0435, 0.045, 0.0465, 0.652,0.0515, 
$0.049, .046, .045, .045,. 046, .047,. 051, .054,, 055,. 057, .044, .046, 
$0.050, 0.051, 0.050, 0.047, 0.044, 0.0435, 0.044, 0.045, 0.04 7, 0.050,0.054 
$,0.056, 0.057, 0.045, 0.047, 0.051, 0.053, 0.051, 0,046, 0.043, 0.042, 0.042 
$,0.0432,0.047,0.050,0.054,0.057,0,056/ 

QATA(CTXCP4(I.J) ,1*1, 15), J*l, 5 )/0.90,0.75, 0.55, 0.38, 0.50,0. 7 2, 0.90 
$,0.94,0.90,0.60,0.60,0.41,0.25,0.15,0.15,0.55,0.40,0.05,-0.10,0.15 
$,. 47, 0,70,0. 85, 0.8S, 0.83, 0.60, 0.37, 0,15, -0.05, -0.1 1,-0. 11, -0.45, 
$-0, 75,-0, 90,-0. 65,-0.25, 0.25, 0,45,0, 53,0. 55, 0.50,0. 30,0. 10, -0.15, 
$-0, 35,-0 ,30,-0. 60,-1 ,20,-1 .30, -1.05, -0.60, 0.0, 0.30, 0.42, 0.40,0.40, 
$0. 20, 0.0, -0.24, -0.45, -0.55, -0.80, -1.46, -1.65, -1 .50,-1 .0,-0. 30,0. 13 
1$, 0.28, 0.39, 0.30,0. 10,-0. 15,-0.35,-0.60/ 

RAD  *  57.29578 
IF  (N0SE.EQ.2)  GO  TO  10 

..................................  OGIVE-CYLINDER 

CALL  TABLE2  12 ,CNA1 , TCNA1 , FMACH , TKACH , 15 , 3 , XLN ,TXLN , 3 , 2 ) 

CALL  TABLE2  (2 ,CNA2 ,TCNA2 , FMACH , TMACH , 1 5 , 3, XLA , TXLA , 3 , 2 ) 

CALL  TABLE2  (2 ,CNAJ ,TCNA2, FMACH , TMACH ,15, 3, 6,0, TXLA, 3, 2) 

CNA  ■  (CNA1+CNA2-CNAJ)*RAD 

CALL  TABLE2  (2,XCP1 , TXCPl ,FMACH, TMACH , 15, 3, XLN ,TXLN , 3 , 2 ) 

CALL  IABLE2  ( 2, ACP2,TXCP2, FMACH ,TMACH , 15, 3 , XLA , TXLA , 3 , 2 ) 

CALL  TABLE2  (2 ,XCP3,TXCP2, FMACH, TMACh, 15, 3, 6,0, TXLA, 3 , 2 ) 

XCP*XLN- ( XCP1+XCP2-XCP3) 

GO  TO  11 

..................................  CONE-CYLINDER 

10  ANG*ATAN  I.5/XLN) 

CALL  7ABLE2  <2,CNA4, TCNA3, FMACH, TMACH, 15,3, ANG, TANG, 4, 2) 

CALL  TABLC2  (2,CNA5, ICNA4, FMACH, TMACH , 1 5 , 3, XLA, TALA, 5, 2) 

CALL  TABLE2  (2,CNA6, TCNA4, FMACH, TNACH ,15, 3,6,0, TALA, 5, 2) 
CNA*(CNA4+CNA5*CN A6 ) 6RAD 

CALL  TABLE2  (2 ,XCP4,TXCP3 ,FNACH , TMACH ,15, 3 , ANG, TANG, 4,2) 

CALL  TABLC2  (2,XCP5,TXCP4 , FMACH , TMACH, 15, 3 , XLA, TALA , 5, 2 > 

CALL  TA0LE2  (2, XCP6, TXCP4 , TMACH, TMACH, 15, 3, 6.0, TALA, 5, 2) 
XCMXLN-(XCP4tXCP5-XCP6)  87 


11  RL1URN 
END 

SUBRUUTINE  BO  AIL  ( FM ACM , DBUD , XLBUU , CN A , XCPOD ) 
COMPUTES  CCNA)  AND  (XCP)  FUR  CONICAL  bOATTAlLS 
REF.  MICOM  REPORT  Nil.  KU-TM-b«-5 
6J.  w.  D.  aAShlNGION 
INPUTS 

FMACH  =  FH EE  STREAM  MACH  NUMBER 

DBUD  *  DIAMETER  RATIO  (BOAT  I A1L/C YL1 NDbK ) 

XLBOD  *  BUA1IA1L  LENGTH/CYLINDEH  DIAMETER 


DBUD 

XLBUU 

OUTPUT 

CNA 

XCPOD 


CNA  *  BOATTA1L  NORMAL  FURCE  CUEFF , ,  bASED  UN  CYLINDER  mKEA-1/KAD 
XCPOD  *  BOATTA1L  CENTER  OF  PRESSURE,  CALIBERS  FROM  LEADING  EDGE 
DIMENSION  TFMACHU  1  )  ,  TaCPLBI  1 1  )  ,  TbETLB  i  1 8)  ,  TCNPf<  (  1  8  ) 

DATA  ITXCPLBII),  1*1, 11)  /. 41, .41, .41, .415, .425, .45, .48, .62, .58, 

♦  . bi5 , . 7 / 

DATA  ITFMACHHJ,  1*1,11)  /  .  0  ,  .  5 , 1  .  ,  1  .  5  ,  2  .  ,  2 . 5 , 3  .  ,  3 . 5 , 4  .  ,  4  .  5  ,  5  .  / 

LATA  (TbETLB(I),  1*1, IB)  /-2  .  , -1 .  5 , -1 .  ,  -  .  7S ,  -  .  5 , -.  2S ,  u . ,  .  25 ,  .  5 , 

4  .7S,1., 1.25,1.5,2. ,3. ,4. ,5, #6./ 

DATA  (TCNPK(l),  1*1, 18)  / 2 . 01 , 2 . 01 , 2 . 03 , 2 . 1 2 , 2 . 52 , 2 . 92 , J. UP , 2 . B 7 , 

♦  2. 35, l.bb, 1.15,1., .93, .83, .69, .57, .5, .43/ 

UETLB*SURT ( ABS I F'M  ACH**2"1 . ) ) /XLBOD 
1FIFMACH.LT.1.)  BETLBs-BETLB 
IFlBETLb.LE.b.)  GO  TO  10 

WRITE  lb, 500)  BET  LB 

500  FORMAT  (1H0, *BETLB*«F10, 2, 5X, •OUTSIDE  TABLE  5.0,  LAST  POINT  USED') 
BETLB*6.0 

10  XMACHaFMACH 
IFIFMACH.LE.5.)  GO  TO  11 
WRITE  lb, 501)  MACH 

501  FORMAT  (1H0, 'FMACH, 'Fil». 2, 5X, 'OUTSIDE  TABLE  5.0,  lAST  POINT  USEu') 
!  XMACH*5.0 

11  call  1NTERP  (2,1 ,CNPR , TCNPR , BETLB, IBETLB, 1 8,2,M1NX,MAXX) 

CALL  INTERP  (  2  » 1 ,  XCPLB,  TXCPLB ,  XMACh ,  TF  *  ACH  ,  1 1 , 2  ,  « I  NX  ,  M  AXX  ) 
CNA*-CNPR4(1.-DB0D**2) 

XCPOD*XCPLB*XLBOD 
I  RETURN 
!  END 

i  SUBROUTINE  CR05FL0 ( FMACH , ALP , RE , D , LUDT , LODN , NOSE , XCG , C* V , CM  V , XCP V ) 

!  CALCULATE  VISCOUS  COMPONENTS  UF  NORMAL  FORCE  AND  PITCHING  MOMENT 
|  USING  ALLEN'S  METHOD  wITH  MODIFIED  CDC  CURVE. 

1  INPUTS I 

I  FMACH  *  FREE  STREAM  MACH  NUMBER. 

ALP  *  ANGLE  OF  ATTACH  -  DEG. 

RE  *  FREE  STREAM  REYNOLDS  NO.  •  1/FT. 

D  a  CYLINDER  DIAMETER  •  FT. 

LOOT  *  TOTAL  BODY  LENGTH  -  CALIBERS. 

LODN  *  NOSE  LENGTH  -  CALIBERS. 

NOSE  ■  1  OGIVE 

*  2  CONE 

XCG  *  CENTER  OF  GRAVITY  -  CALIBERS  FROM  NOSE. 

OUTPUTS} 

CNV  ■  VISCOUS  COMPONENT  OF  NORMAL  FORCE 

CNV  *  VISCOUS  COMPONENT  OF  PITCHING  MOMENT,  ABOUT  XCG. 

XCPV  a  VISCOUS  COMPONENT  CENTER  OF  PRESSURE  -  CAL.  FROM  NOSE. 

NOTES} 

1  CALIBER  IS  l  CYLINDER  DIAMETER  (0). 

CNV  AND  CNV  ARC  BASED  ON  CYLINDER  CROSS  SECTIONAL  AREA 
AND  DIAMETER. 

DIMENSION  TL0D(7 ) , TETAOC 7 ) , TMN (22 ) , TCDC ( 22 ) ,TMN2(6), TCDC2 (6) 

REAL  MN, LOON, LOO A, LOOT 

OATS  TLOO  /O., 5., 10. , 15. ,20., 30., 40./ 
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DATA  IfcTAU  /.Di,  ,bi,  .69,  .  74,  ,7b,  ,«0,  .*2/ 

DATA  TMW  /,  0,  .  1  ,  .2,  .  j,  .  4,  ,5,  .6,  ,7  ,  ,6  ,  ,H5  ,  .9,  .94  , 

^  liUfltlrl|2rla4|l|brlftif2«Of2a4fitbf3«0/ 

DATA  I CUC  /1.2,l.l*,l.ltf,1.2,1.2J,l.2H,l.J7, 

♦  1 .52,1 .45, 1.42, 1.64, 2.08,2.0, 

4  l.b7,l.UU,l.d2,l.B5, 1.85, 1.85, 

4  1.8,1.69,1,50, 1.49,1, 43, 1.40, 1.3b, 1.34, 1.34/ 

DATA  1MN2  /.0, .1 , .2, .3, .4, .5/ 

DATA  rCDC'2  /. 30,. 3  7,. 4b,, 72,. 9t>,1.28/ 

DEG*57 . 2957  8 
PI*3. 14159 
SHALt*l ,e-30 

IF  (NUSfc.EU.l)  GO  TO  16 

CONE  PLANFuKrt  AREA  ANU  CLMHUIU 
APl»(LUDN/2.)*4./Pl 
XB1».6667*L0UN 
GU  TU  15 

0G1V£  PLANFURM  AREA  ANU  CENTROID 
16  bK0D3L00N**2  4  .25 

XB1*,65*L0DN 

APl*-2.*LODN*(bROD-.5)  4  L0DN4SQKT(BR0D**2“LUDn*42 ) 

4  4BR0D**2*ASIN( LODN/bHUD J 
AP1»AP1*4./PI 

ADD  CVLlNDEKlCAL  AF  TfcKBOD T 
15  LOUA*LODT  -  LODN 

AP2*U0DA*4,/P1 
XB2».5*LODA 
APUARsAPl  4  AP 2 

XCPV«(XB1*AP1  4  (LOON  4  XB274AP2J/APGAR 

MN*FRACH*S1MALP/DEG) 

RfcN*KE*D*SlN( ALP/DEG 1 

CALL  1NTERP  (2, 1  .ETAO.IETAO,  LUDT  ,TL0D  ,7,2,*IN,MAX) 

RECR1T*1.£05 

ET A*ET AO 

IF  (REN.Lfc.RfcCRlTl  GU  TO  10 
ir  CNN. GE. 0.5)  GO  TO  10 

CALL  1NTERP  (2,l,CDC,rCOC2,NN,TMN2,6,2,PlN,"AX) 

GO  TO  11 

10  CALL  1NTERP  ( 2 , 1 ,CDC,TCDC ,MN ,TNN , 22 , 2 , MIN , MAX) 

11  CNV*CUC*ETA*AP0AR*S1N( ALP/DEG) **2 
CMV*CNV*(XCG«XCPV) 

RETURN 

END 

SUBROUTINE  CROSFL1  ( FM ACH , ALP , LT , NOSE , XNOSE , XTOT , XCG , CN ^ , CM V , XCP V ) 
CALCULATE  VISCOUS  COMPONENTS  OF  NORMAL  FORCE  AND  PITCHING  MOMENT 
USING  CROSS  FLO*  DRAG  METHOD  OF  PERKINS  AND  JORGENSEN. 

INPUTS! 

FMACH  «  FREE  STREAM  MACH  NO. 

ALP  *  ANGLE  OF  ATTACK  -  DEG. 

LT  «  1  LAMINAR  BOUNDARY  LATER. 

*  2  TURBULENT  BOUNDARY  LATER. 

NOSE  ■  1  OGIVE  NOSE. 

■  2  CONE  NOSE. 

XNOSE  ■  NOSE  LENGTH  •  CALIBERS 

XTOT  ■  TOTAL  BOOT  LENGTH  •  CALIBERS 

XCG  *  CENTER  OF  GRAVITT  -  CALIBERS  FROM  NOBE 

OUTPUTS! 

CNV  ■  VISCOUS  COMPONENT  OF  NORMAL  FORCE 
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*  VISCOUS  COMPONENT  UK  PITCHING  MOMENT 


r 

CmV 

notes : 

l  CALI DEM  IS  1  CTLINUtP  DlAMETEH  (D). 

CN*  AND  CM V  ARE  bASED  UN  CYLINDER  CROSS  SECTIONAL  AREA 
And  DIAMETER. 

DIMENSION  TXOLmI  21 )  , TCOkL(21  J  #  TCDKT(  21  ),TMC(21),TCDCP121) 
DIMENSION  TdUXNU  1  )  ,  TLMUl  1 1  )  ,  '1LMC(  1 1  ) 

HEAL  Ln,LmUAN,MC 

DATA  IaOL*  /.0,.1..2,.i,.4,.S,.b».7,.B,.9,1.0,l.i,1.2, 

♦  1 • 3 , 1 . A  *  I • 5 , 1 • b , 1 . 1 , 1 • » , 1 • 9 , 2 »D/ 

DATA  TCUnL  /  .  0 ,  ,  025  *  .  Ob  5 ,  .  1  H5 ,  .  35 ,  .  b3  ,  .  9  3  , 1 . 2 1> ,  1  .  38  ,  1 . 49 , 

♦  1.48, 1,40, 1.30,1, 21,1.14, 1.08, 1,04, 1.00, .97, .94, .92/ 

OATA  TCDHT  / , 0 , . U 1 , . Ob , . 1 0 , . 25 , . SO , . HO , 1 . 0 , 1 .  1 , 1  .  I , 

♦  1.03, ,9b, .90, .85, .82, .80, .80, .80, .80, .80, .80/ 

DATA  INC  /.O, . 1 , .2, .3# .4, .5, .0, ,7, .8, .9, 1 .0, 1 . 1 , 1 .2, 

♦  1 . 3 ( 1 .4, I • 5 , 1 »b , 1 . 7 , 1 , 8 , 1 . 9 , 2.0/ 

DATA  TCDCP  / 1 . 2 , 1 . 1 8 , I • l b , 1 . 2 , I . 2b , 1 . 3b , 1 . 53 , 1 . 7 3 , 1 .  82  , 

♦  X.8'2, 1.79,1.  74,  l.b7,l  .bO,  1.53.1 ,4b,  1.40, 1. 37,1.33, 1  .  30,  l 
DATA  TbuXN  /.Q,.l,.2,,3,.4,.5,.b,.7,.B,.9,1.0/ 

DATA  TLMO  /I ,0 , 1 .09, 1 . 22, 1 , 39, l . 57 , l . 7b , 1 . 98 , 2. 20, 2 . 42 , 2. o5 
DATA  TLMC  /I ,0, 1 .21 , 1 .42/ 1 ,b3, 1 . 84 , 2 .Ob, 2. 27 , 2.4  1 ,  2  .b  1 , 2.89 

SMALL*1,E»30 
DEG*57. 2957b 
Pl*3. 14159 

LMUXNsl . 0 

IF  IFMACH.LE. 1.0)  GO  TO  11 
DU AN  SSURT l I MACH**2-1.)/XN0SE 
IF  (NUSE.EU.2)  GO  TO  10 

'CALL  IN1EKP  (2,1,LM0XN, TLMU,  BOXN  ,1'BUXK  ,  11 ,2, MIN, MAX) 

GU  TO  11 

10  CALL  INTtRP  (2,1,LM0XN,1LMC,»UXN,1B0XN,11,2,MIN,4AA) 

11  lm»lmoxn*xnose 

! . .  CALCULATE  CDCP 

MCSFMACH*S1N( AbP/DEG) 

|CALL  INTERP  (2,1,CDCP,TCDCP,mC,TMC,20,2,MIN,MAX) 


■S  U  m  N  a  0  , 

SUMM*0. 

FNlaO. 

FM1»D, 

X1*0. 

1*2 

56  X2*X 1  4  ,1*LM 

ir  (X2.LE.XT0T )  GO  TO  54 

X2*XT0T 

X0LM*X2/LM 

IF  (LT.EU.2)  GO  TO  55 

CALL  INTERP  ( 2 # 1 »CDR2 » TCDRL » XQLM , TXOLM , 21 , 2 , MIh , M AX ) 

GO  TO  51 

55  CALL  INTERP  (2 , 1 , CUR2 ,TCDRT, XOLM ,TXULM , 21 , 2 , MIN , MAX ) 
GO  TO  51 

54  IF  (LT.EQ.2)  GO  TO  50 

CDR2*TCDRL(1) 

GO  TO  51 

50  CDR2*TCDRT 11) 

51  M2* ,5 

IF(X2.GE.XN05E)  GO  TO  52 

ir  (NU5E.EU,2)  GO  TO  52 

•U»R*.2S  4  XN0SE442  90 


28/ 

2.8«/ 

3.U9/ 


H  2*SUK  T  (  H 1  GR ♦  * 2- (  XNOSfc-Xi  )  *  *2  )  •  ( BI  GR- .  5 ) 

GO  TU  Si 

52  K2ax2*U./XNOSE')/2, 

53  FN2aCDK2*R2 
FM2»Fn2*(XCG-X2) 

0N»(FN2  4  FNn*lX2-Xl}/2. 

0Ma(FM2  4  FM1)4(X2-X1  )/2. 

SUNNaSUMN  4  UN 
SUMMasUMM  4  Om 

Xl»X2 
fN  t*FN2 
FMlaF*2 

1*1  4  l 

IF  (1.GT.21)  1  =  21 
IF  (X2.LT.XTOn  GO  10  56 

CONST* (8. /PI )*CDCP*S1N( ALP/DEG)**2 

CNV*CUNST»5UMN 

CMV*CUN5T*SUMM 

XCPVaXCG  -  CMV/(CNV  4  SMALL ) 

RETURN 

ENO 

SUBROUTINE  ARATIXM, CCR  ,  XX  1  , XALB, DDb ,  D  ,  A 1  A  ) 

IHIS  SUBROUTINE  calculates  THE  SHADED  AREA  RATIO 

FOR  HODIFING  CARRYOVER  »RE»  FINS  ARE  CLOSE  TO  BASE  OF  BODY 

FUR  CYLINDER  OR  BOA1TAIL 

XH  a  NACH  NUMBER 

CCR  a  FIN  ROOI  CHORD 

XXI  a  DISTANCE  FROM  BASE  OF  MISSILE  TO  FIN  TRAILING  EDGE 
XXLB  a  LENGTH  UE  BOAT1A1L 

ODB  a  DIAMETER  OF  BASE 

0  a  DIAMETER  OF  CYLINDER 

A 1 A  a  SHADED  AREA  KAliU 

PHUGRAM  MODIFIED  19  DEC, 71)  10  INCLUDE  FINS  EXTENDED  AFT  BASE ,  IAS 
.  A1A*1.0 

,  IF(XM.LE.l.O)  GU  TO  7 
CR*CCR/D 
X1*XXI/D 
XLBaXXLB/D 
DBaDDB/D 
AaCR 

IFtXl.LT.  0.0}  CRaCR4Xl 
IFU1.LT.  0,0)  XIaO.O 

XLCaCR  4  XI  -  XLB 
B*TAaS0RT(XM*XM  -  1.0) 

TANMUal. O/BETA 
TANTHaO.O 

IFIXL8.NE.0.0)  TANTH*,5*(1.0-DB)/XLB 
TA»TM»(TANTH  ♦  TANMU}/( 1.0-TANTH*TANMU) 

C1*S0RT( 1 ,0  ♦  XLC*XLC) 

C2*SQRT (1.0  ♦  ( (XLC  4  XLB)/ l .5* ( 1 .0  ♦  DB)))**2) 

C3«S0RT(1,0  ♦  (tANTH  4  Xl * l 1 , 0  t  TANTH**2)/DB)**2) 

IFCCXM.GE.  1.0).  AND,  IXM.LE.CD)  iai 
IF(  (XM.CT.Cn. AND, (XM.LE.C2))  1*2 
2F( (XM.GT.C2 ) .AND. (XM.LE. 1 .0E42S) )  1*3 
2F( (XM.GT.C3) . AND, ( XM.LE, 1 .0E425) )  Jal 
IF(CXM.CE.1.0),AND.(XM.Lb.C3))  Ja2 
60  TO  (1,2,3),  1 
t  8*0, 
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H*BETA 

C*0. 

H  1*0  , 

GO  TO  (4,5),  J 

2  b=l .O-XLC/BEIA 
HsXLC 

C*0. 

ril*b/(UMn  ♦  TAnMU) 

GO  TO  (4,5),  J 

3  bsI.O-XLC/bEIA 
H*XLC 

C*.5*( 1  ,0  +  0bJ-UtiC4XLtt)/tJt'lA 
H 1 sXLb 

GO  Tu  (4,5),  J 

4  fcal  ,0-2.0*(XLB-Xl  )4TAMH 
FsDb-Xl*(TANlM-rANTH) 

H2*X1 
GO  TO  b 

5  E*1.0-2,0»lXLb-Xl)*l'ANrH 
F*0, 

H2*E/(TANTH-H4NTM) 

GO  TO  6 

6  A1*,5*(2.0»X0C+(1 ,0*DB)»XLB-( 1 .  0  +  B  )  *H- (  b  +  C  )  *H  1  -  (  E  *0e )  *  A 1  ■*•  (  h  +  F  )  *H2  ) 
A1ASA1/A 

7  RETURN 
END 

SUBROUTINE  KCAOX  ( UU , UKECL , HOD , f UU , LH , Lb , EUFF J 

IHPKUVED  VEKSIUN  OF  FFKEAD 
B i  J.M.  CHAF1N 

DIMENSION  FURM(b) , IBUF(132) ,JBUF(i3) , HOb ( 1 ) , FLO ( 1 ) 

1060  FORMAT  (132A1) 

1260  FORMAT ( ’ (A* ,12, 1 ) 1 ) 

1270  FORMATt ’ C ,12, 'XA • ,12, ')' ) 

1260  FORMATt * (' ,12, *A4,A» ,12, ' ) ' ) 

1290  FORMAT ( ' l • ,12, 'X* ,12,  'A4,A* ,12, • ) • ) 

1360  FORMAT ( ' ( • , 12 , • XF • , 1 2 , ' , 1 ) '  ) 

1400  FORMATt ’ tF* ,12, ' .1) ' ) 

1460  FORMATt ' t • ,12, 'XI  * ,12, ) 

.  1500  lORMAlt • tl • ,12, » ) ' ) 


PURPOSE 


REAOS  ONE  CARO  IN  FREE  FIELD  FORMAT  (UP  TO  132  CHARS) 

ANY  COMMENT  CAROS  ENCOUNTERED  IN  THE  DATA  DECK  ARE 
SNIPPED  AND  NOT  INTERPRETED 

ARGUMENTS  IslNPUT 

— — —  u*ourpur 

LU  —  1  —  LOGICAL  UNIT  NUMBER  OF  INPUT  DEVICE 
LRECL—  1  —  LOGICAL  RECORD  LENGTH  IN  CHARACTERS  U32  MAX) 

HOL  —  0  —  MILL  CONTAIN  HOLLERITHS  FOUND  ON  RECORD  t33  MAX) 

FLO  —  0  —  NILL  CONTAIN  FLOATS  FOUND  ON  RECORD 

LH  ••  0  —  NUMBER  OF  HOLLERITHS  ON  RECORD 

LF  —  0  —  NUMBER  OF  FLOATS  ON  RECORD 

COrF  —0—0  RETURNED  IF  NO  END  OF  FILE  ENCOUNTERED, 

corr  —0—1  RETURNEO  ir  END  OF  FILE  ENCOUNTERED, 

DEFINITIONS  —  LIMITATIONS 
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A  ♦TRACAK*  l  TRANSPARENT  CHARACTER)  IS  A  bLANK,  COMMA  UK  a  SIGN 
THE  FOLLOWING  AKE  DELIMITERS 

-  BEGIN  UF  CARD 

•  END  UF  CARO 

*  A  SIRING  CONSISTING  OF  ONE  OK  MORE  TRACARS 

AN  ITEM  CONSISTS  OF  ONE  OK  MOKE  NUN-TRACAKS  PRECEDED  AND  FOLLUmED  ol 

DELIMITER. 

EACH  ITEM  RILL  b£  INTERPRETED  AS  EITHER  INTEGER,  FLOAT,  OK  HOLLERITH 
ANT  ITEM  STARTING  WITH  +  -  .  UK  DIGIT  KILL  BE  INTERPRETED  AS 

-  INTEGER  IF  IT  CONTAINS  NO  . 

-  FLOAT  IF  IT  CONTAINS  ONE  . 

-  HOLLERITH  IF  IT  CONTAINS  MORE  THAN  UNE  . 

ALL  OTHER  ITEMS  «1LL  HE  INTERPRETED  AS  HOLLERITHS 


INITIALIZE  COUNTS  AND  READ  CARD 


£OFF*0 

1020  LHSO 
LFao 
LI»0 

1021  REA D(LU, 1060, £ND*999)  U BUF ( I ) , lai , LRECL  ) 
IFdbUFd)  .  £Q  *  •  C  *  .AND.lbuFl2I.EQ.  1  ')  GO  TO  1021 
ENCODE IJBUF, 1060)  ( IBUF ( I ) , lal , LRECL ) 

12*0 


HMD  II  a  FIRST  COLUMN  OF  FIELD 


1100  Ilal2 
1120  11*1 1*1 
IF  (ll.GT.LKECL)  RETURN 
lXalbUF(ll) 

IF  (IX.EU. '  • .UR. IX.EU. .OH. IX.EU. •*' )  GU  TO  1120 


FIND  12  a  LAST  COLUMN  OF  FIELD 


12*11 

1140  12*12+1 

IF  (12. GT. LRECL)  GO  TO  1160 
IX*1BUF(I2) 

IF  (IX. NE.*  ’ .AND.IX.NE. ' , * .AND.1X.NE, '*' )  GO  TO  1140 
12*12*1 
GO  TO  11«0 
1160  I2*LRECL 


IN  •  FIELO  mIOTH 


U«0  NX*I  1*1 
IN*. 2-11+1 
00  1190  1*1,6 
1190  rORN(I)  «'  * 

IX*1BUF(I1 ) 

IFdX.EO.'  +  '.OK.IX.CU. '-’.OR.IX.EO. GO  TO  1200 
ir  (IX.CC.'O'.AND.IX.LE. '9')  GO  TO  1200 
GO  TO  1240 


COUNT  THE  DOTS 


1200  NDOTS*0 

00  1220  1*11,12 

ir  dbUFCD.EO.'.  ' )  N DOTS* N DOTS +  1 
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1 

122U  CuMlNUE 

lF(NOOTS-l)  1440, 1140, 1240 


HGLLbhlTH 


1240  CONTINUE 
nF*1*/4 

NP*NAX0(1 ,MUD(1*,4)  ) 

IF(N1 .EU.O.AND.NX.EO.O)  ENCODE (FORM, 1260) 

IF(NF.EU.O.AND.nx.GT.O)  encode i FOR M, 12 ?0) 

IF(NF.GT.O.AND.wX.EQ.O)  ENCODE  (FORM,  12  WO) 

IF(NF.GI.O.AND.NX.GT.O)  ENCODE  (FORM,  12 'JO) 

1 F ( NOD ( 1 « , 4 ) . Gi . 0 )  NF*NF+1 

DrtlsLn+1 

LHN*DM+NF 

DECODE  (JBU1  ,FURM)  (HOLU  )  ,  1*LH1  ,LHN) 

LH*Lm*NF 

GO  XU  1100 


NP 

NX,NP 
NF  ,  N P 
NX , NF , NP 


FDUAT 


1 J40  CONTINUE 

IF(NX.EO.O)  ENCODE(FORM,  1400)  1*. 
IF(NX.GT.O)  ENCUDE(FUKN , UbO)  NX,I« 
LF*LF+1 

DECODE  (JBUF.FORM JFLO(LF) 

GO  TU  1100 


INTEGER 


1440  CONTINUE 

IF(NX.fcO.O)  ENCODE (FORM, 1500}  1* 
IF(NX.GT.O)  ENCODE  (FORM,  14  00)  NX.1* 
LF*LH-1 

DECODE  (JbUF,FORM)  XNT 
FL01DF)*FLUAT(1NI) 

GO  TO  1100 


IF  END  OF  FIDE  ENCOUNTERED,  SET  EUFFsl 


*99  EUFF*1 

RETURN 

END 

SUBROUTINE  ROLL  ( BO , CR , TR , KBODY , DR EF , AREF , CLP , CLU ) 

SUBROUTINE  FOR  CALCULATING  ROLL  DAMPING  AND  ROLL  MOMENT  EFFECTIVENESS. 
THEORY  FROM  nACA  REPORT  1088, 

ALL  GEOMETRIC  INPUTS  MUST  HAVE  SAME  SYSTEM  or  UNITS. 

I.E.,  FT,,  IN,,  CALIBERS,  ETC, 

BO  «  FIN  EXPOSED  SEMI-SPAN 
CR  ■  FIN  ROOT  CHORD 
TR  *  FIN  TAPER  RATIO 
RBODY*  RADIUS  OF  BODY  AT  FIN 
DREF*  REFERENCE  LENGTH 
AREF*  REFERENCE  AREA 

OUTPUT 

CLP  «  ROLL  DAMPING  COEFFICIENT  (BASED  UN  DREF  f  AREF) 

CLD  *  ROLLING  MOMENT  EFFECTIVENESS  COEFF  (BASED  UN  DREF  $  AREF) 

(FOUR  FINS  DEFLECTED) 
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UREFskEF EKEnCE  DIAMETEk 

DIMENSION  ILAMDAI  n  )  ,  TCLPAkl  1  1  )  ,  ICLOARU  1  ) 

DATA  TLAMDA/,0, .1 , .2, .3,.4,.5,.t>,.7,.«,.9,1.0/ 

DAIA  TCLPAR/.16, .159, .161 , .164, ,loU,.1475, .120, .u»b, .U4J, .01 J, .u/ 
DATA  XCLUAK/.252,.254,.25o,.24tt,.23Q,,202,.15b,.100,.044,.014,,0/ 
CLPau.O 
CLDau.O 

DB0C*»2.4RBUDY 
IF ( du . EU .0,0}  KETUKN 
SPAn«2.»B0  ♦  D6UD i 
XLAMLAsDBODY/SPAN 

SaCk*  (  BQ+DbOQ!  /  2  ,  )  4  (  1  .♦IK+ll  .-Tk)*DBUD)r/(2.*B<J)  ) 

AR*SPAN**2/S 

CALL  I^TEkP  (2,l,CLPAR,TCLPAk,XLAMDA,TLAMDA,ll,3,*UN,M4X) 

CALl  I  n  T  E  k  P  (2, 1 ,CLOAk,TCLDAR,XLAMUA,lLAMDA, 11 , J,M1N,MXX) 

CLP*-CLPAk*AK 

CLO«CLUAK»AK 

CLPsCLP»CSPAN/UR£f )**2*S/AKEf 

CLUaCLD»lSPAN/UKEF)*S/AREF 

RETORN 

END 

SUBROUTINE  FINGEOM  (CR,TR,BQ,AR£F,SWTE,CT,AR,S*l£,CBAH, 

♦  T BA R,XM, CENT, AFAR, TANL,S«HA) 

FIN  PANEL  GEOMETRY. 

ALL  GEOMETRIC  INPUTS  MUST  HAVE  SAME  SISTER  UF  INPUTS. 

I.E.,  FT.,  IN.,  CALIBERS,  ETC. 

OUTPUTS  MILL  HAVE  CONSlSlANi  UNITS  wlTH  INPUTS. 

INPUTS: 

CR  *  ROOT  CHORD 

IR  *  TAPER  RATIO 

BU  a  EXPOSED  SEMI-SPAN 

AREF  «  REFERENCE  AREA 

SATE  *  TRAILING  EDGE  SmEEPBACK  ANGLE  •  DEG. 

OUTPUTS: 

CT  *  TIP  CHORD 

AR  a  ASPECT  RATIO  1  1  FIN  PANEL  }. 

(  SAME  FOR  2  FINS  JOINED,  WITHOUT  BODY), 

SwLC  a  LEADING  EDGE  SwEEPBACK  ANGLE  -  DEG. 

CBAR  a  MEAN  AERODYNAMIC  CHORD  /  ROOT  CHORO, 

YBAR  a  5PANW1SE  LOCATION  OF  M.A.C,  /  EXPOSED  SEMI-SPAN. 

XM  a  DISTANCE  FROM  L.E.  CR  TO  L.E.  MAC. 

CENT  *  DISTANCE  FROM  L.E,  CR  TO  CENTROID. 

AFAR  a  AREA  OF  1  FIN  /  AREF. 

TANL  a  TANGENT  OF  LEADING  EDGE  SNEEPBACK  ANGEL 

S«HA  a  MID  CHORO  SwEEPBACK  ANGLE  •  DEG, 

SMALLal.E-iO 

RA0«S7.29S7« 

ARa4.*B0/(CR*ll.+TR) ♦SMALL) 

CT»CR*TR 

AFAR«CR*B04 ( 1 . ♦TR ) /AREF/2 . 

TANLa(CR-CT)/(BUtSMALL)tTAN (SwTE/RAD) 

SBl,E»l  ATAN  (TANL))*57.3 
TARHaTANL*. 5* (CR-CT)/ (BO+SMwLL) 

SNNAa t ATAN (TANH) ) *57 , 3 

CBARal .0-( 1,0-(2,0*TR**2)/(1,0+TR) )/l,0 
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YbAR«U(t.042.0*TR)/(l.G+TR)  J/3.0J 
X**loAR*I ANL*BU/ (CR+SMALL) 

CENT*X**,5*CbAR 
RETURN 
END 

SJbRuuTUt  SI AN62 (HEIGHT, VSOUNU,KHU, PRESS, TEMP, I UNITS) 

THIS  SUBROUTINE  IS  A  U.  6.  STANDARD  MODEL  ATMOSPHERE,  1462 
PREPARED  BY  AMSMl-HST,  21  JULY  19b7 

THIS  ROUllNE  IS  VALID  BELUn  700,000  METERS 
EXTRAPOLATION  IS  PROVIDED  ABUVt  700,000  METERS 

1  ntlGHl  ALTITUDE  M 

l  V SOUND  VELOCITY  OF  SOUND  M/SEC 

i  KHu  DENSITY  KG/MJ 

4  PRESS  PRESSURE  KG/M  2 

5  TEMP  TEMPERATURE  K 

e  UNITS  >1  FUR  METRIC 

*2  FOR  ENGLISH 

DIMENSION  ZZ(li),PB(l33'RLll3>»ZB(133,lb(lJ)  ,XXX(liJ 
DIMENSION  HH(8) ,PBH(8),RLH(8),HBH(B) ,TBH(B) 

DATA  HH/  11.0,  20.0,  32.0,  47.0,  52.0,  bl.O,  79.0,  86.743/ 

DAIA  PBH/  10331.907b,  2307.7398,  558.2615,  88.50965,  11.30876, 

1  0.016163,  1.85682,  .105812/ 

uAIA  RLH/  >6.5,  0.0,  1.0,  2.6,  0.0,  -2.0,  -4.0,  0.0/ 

DATA  HBH/  0.0,  11.0,  20.0,  32.0,  47.0,  52.0,  61.0,  79. J/ 

DATA  TBH/  268.15,  216.65,  216.65,  228.65,  270.65,  27u.65,  252.65, 

1  160.65/ 

DAIA  11/  100. , 110. , 120. , 150. , 160., 170. ,19u., 230. ,300. ,400. ,500., 

1  600, ,700./ 

DATA  PB/  .016762,  .0030067,  .7499E-3,  .2571E-3,  ,51oE-4, 

1  .3767E-4,  . 28476E-4,  .1/1B4E-4,  .709767E-5,  .1921E-5,  .41U97E-6, 

2  .1117E-6,  .3518E-7/ 

DATA  RL/  3.0,  5.0,  10.0,  20.0,  15.0,  10.0,  7.0,  5.0,  4.0,  3.3, 

1  2.6,  1.7,  1.1/ 

DATA  ZB/  90.0,  100.0,  110.0,  120.0,  150,0,  160.0,  1 70.0,  190. 0, 

1  230.0,  300.0,  400.0,  500.0,  600,0/ 

DATA  IB/  180.65,  210.65,  260.65,  360.65,  960.65,  1110.65,  1210.65, 
j  1  1350.65,  1550,65,  1830.65,  2160.65,  2420.65,  2590.65/ 

I  DATA  XXX/  33.166119,  33.064396,  32.96512,  32.784318,  32.530427, 
i  1  32.455229,  32.309443,  32.016142,  31.492584,  30.704189, 

1  29.808321,  28,937691,  28.123252/ 
i  DATA  AA,BB,CC,DD,£E/ 

1  1 . 5731 262E-7 , 2 . 4656553E-14, 3 . 8667054E-21 ,6,0621 354E-2 8 , 

2  9, 501364 VE-3S/ 

«MMII*MMM**M*tl****MMMM*  0  -  90,000, 


ir  (height.lt. o.ooon  height«o.o 
H  m  height 

iPdUNiTS  .ME.  1)  H  *  H  ♦  .3049 
Z»H*,001 

I rcz  .GT,  90.03  GO  TO  140 

H  •  (H  •  AA*H**2  ♦  BB*H**3  •  CC*H**4  ♦  DD*H**5  -  CE*H«*6)  *  ,001 
00  S50  1*1 , t 

irCH  •  HH( 1 3 )  551,551,550 
950  CONTINUE 
GO  TO  140 

551  1F((1.NE.2) . AND. (1 ,N£.5) .AND. (1 ,NC. 8) )  GO  TO  120 
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FT 

FT/SbC 
SLUG/f TJ 
LBS/FT2 
H 


ho 

PHtSS*FeH(i)*£XP(-34.1b«79*(H-HBH(l))/ibH(l) ) 

GU  TU  130 

120  I  fc.**P*TBH  H  )  4  (H-HBH(i))*RLH(I) 

EP*34.1b479/RLH(I) 

PRESS*Pd«  U  )  *  ( TBM  ( 1 ) /TEMP )  **EP 
130  Kht  «  PK£SS/(29.2b945*T£MP) 

V60UNL.  s  20.04b333*TEMP**0.5 
GO  TO  4bL 

»«*t*«*»«««***»**«*******«4*******  90,000.  •  700,000 

140  Zl  *  HEIGHT 

IFUl  ,0T.  17 0000. 0)  oo  TU  142 
2  BETWEEN  90000.  *  AND  170000.  M 

XN  *  U  ( ( (0. 141HB509E-2 7*Z1  -0 . 1 1 1  450  34  1 1-2 1  )  *2  l 

1  4U. 35920l41bt-lb)  *21  -0 . bUbb52 1 JL- 1 1 ) *Z 1 

2  t0.bbb21Slb Jt-b)  4  21  -0 . 2  1  5  J0035bE- 1 )  *2 1 

3  45.76957191E42 
GU  TO  150 

2  GKtATER  THAN  1/UOUU.  H 

142  XN  *  (l(((1.10O55*E-33*Zl  -3.4291b2fc-27) 

1  *Zl43.0159bb9£-21 )*Z1  -2 . 0402054E-1 5 ) 

l  *21  ♦5.b422141E-10)*21  -1 .0700489E-4) 

3  *21  43 • 5629995E* 1 

150  00  500  1*1,13 

IF ( 2*22 (1)1  450,450,500 
500  CONTINUE 
*NlTb  lb, 2901 

290  FUHNAT130H0  ALT  ,GT.  700000.0  M..;STANo2) 

1*13 

450  TMsTBlI)  4  (Z-Ztt(i J )*HL(1  J 

fEHP*TM*AM/24.9b4 

LP*XXX(I  )/RL(I  ) 

PKESS*PB(I)*(TB(I)/TM)**EP 
RHo  *  PRESS/ (29.26945*TM) 

V SOUND  u  (401.904b?  *  TM)  **  0,5 

460  IF ((UNITS  .LO.  li  GU  TU  999 

tfSUUNU  *  VSOUNU  *  3,2009333 

RHU  *  KHU  *  .  1940  32E-2 

PRESS  *  PRESS  *  .20402 

TEMP  «  TEMP  *  1,0 

999  RETURN 

END 

SUBROUTINE  BRIT  (XM,AR,TR, SHEEP, CLA,XCPCB) 

TAKEN  FROM  THE  BRITISH  DATA  SHEETS  S. 01. 03. 06 

S. 01. 03. 05 
S. 01. 03. 04 
S. 01. 03. 03 
S. 00. 01. 02 

XN  *  FREE  STREAM  MACH  NO. 

AR  •  ASPECT  RATIO  (T*0  NlMGS  JOINED  WITHOUT  BODY) 

TR  «  TAPER  RATIO 

SHEEP  *  WING  MID-CHORD  SwEEPBACK  ANGLE  ( UEG j 
CLA  •  LIFT  COEFF.  SLOPE,  BASED  ON  EXPOSED  WING  AREA  1/KAO 

XCPCB  *  CENTER  OF  PRESSURE  (PERCENT  MEAN  GEOMETEHIC  CHORD), 

MEASURED  FROM  LEADING  EDGE  MEAN  GEOMETERIC  CHORD. 
DIMENSION  X11(15),YTI(7),2T1(4),XT2(«),YT2(5),ZT2(4),YT4(7) 
DINCNSIUN  BARI (15 , 7,4),  BAR2(« ,5,4) ,BAR3( 1 5,7 ,4) , BAR4( 0, 7 , 4) 
DATA  XT 1/0. , , 5, 1 . , 1 , 5, 2. , 2 . 5, 3  # , 3. 5, 4 , , 4 . 5, S . , 5 . 5 , 6 . , 6 . 5 , 7 . / 
DATA  YT1/0.,1.,2.,3,,4.,S.,6./ 

DATA  211/0. , .25, .5, 1./ 

DATA  XT2/0. , l.,2.,3.,4.,5.,6.,7./ 
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•  7  9  §  0  ^tf#*77#«75#*73#*70#,6W 
52,  1.57,1.53,  1.4b,  1.34, 1.25, 1.1 6, 1.07,1, 
o,  .  7  35,.b95». bb,l. 5 7, 1.53 #1.42, 1.31,1.20 ,1.11 ,1.0. 

l#7o, #74, #70,. bob, .03, 1.2b, 1.19, 1#12, 1*06,1*00, .9! 
i.  »  74#  #  7b#  .  b  7  #  #  b  4 ,  ■  b  1  #  1  .Ob,  1  #00  #  #  9b  #  #917#  #  89b  #  •  4  b  #  #  i 
j>b^,.bbS--K4--h'  <ia  .  .  a»i .  _  i-L  »o  .  mi.  lx.. 

4.b2,.bO,.57,.55,,79,(7b,.73,.71,.70,.b8,.b7,.bb 

<.  m  «>>/ 


01 

.57 

,1.52 

,1.45 

7  1  f 

34, 

1. 

22 

>1.12 

#  1*04#  # 

97,. 

90 

9  • 

9. 

o4. 

1.57 

# 

1.51 

9 

1.42,1 

.30, 

1.18, 

1.08 

,1. 

00, . 

94 

9  • 

1 . 

05, 

.  b3 , 

i 

.57, 

1 

.48 

,1. 

3b 

9  1 

22,1 

.10, 

1.02, .95, 

.8 

2. 

03, 

•  61# 

i 

.27, 

1 

.17 

9  1  8 

10 

9  1 

025, 

.9b, 

.91 

,.8b 

9  • 

81 

3. 

bO, 

•  5b  # 

i 

.Ob, 

• 

579, 

.94 

9  • 

89 

•  84, 

.81, 

.78 

7  • 

9  • 

71 

l)  A  I A  b  AR2  / 

A 

1.347,  .327,.  314,. 
2.31b#. 313#. bOb#. 

3. 410. . 395. .383.. 

4.202. . 265. .267.. 

5.242. . 242. .282.. 

6.365. . 345. .325. . 
7. 24b, .177, .217,. 

8.250. . 250. .248.. 

9.287. . 275. .270.. 

1.245. . 245. .085.. 


Il,^7,l, 56 ,1.85, 1.82, 1.64, 1.41, 1*21,1. (76,  .94,  .845,  .7  7 , 
8.558, 1.57, 1.50, 1.83, 1.72, 1.585, 1.41, 1.215, 1.065,. 955, 

1.66..  61.. 562. 1.57. 1.1 3. 1.46. 1.46. 1.40. 1.305. 1.22. 1.11 

2. 804..  71. .67. .62.. 572. 1.26. 1.18. 1.10. 1.23. 1.275. 1.20, 
31 ,015,. 94, .85, .772,, 71,. 65,. 60, 1.05, 1.05, 1.05, 1.05,1. 


4 1  .  j  ,  i^0Haa9lS^#d^ia«dJa#^42j 

5*9ia  i9i!)f « 9 1 1  »b75  j  ■  H  4 1  iblb  a  «  7 1 3 ,  #73#  ,  b  7 ,  #79#  •  d  b  #  #7(4#  #  7(7#  #7  79#  #79# 
b«7ibi ,(u# 1 47  # .»3, ■ B 7  # , 70b# ,7Sb# #73#  #704/ 

041*  iUbAH3(l»J,K),l*l,lS),J»l,7j,K«3,4J  / 

A  1.58, 1.08, 1.92, 1.80, 1.59, 1.3b, 1.1 7,1. 03,. 92, .04, 

1./5,. 09, ,b35,« 59, .55, 1.50, 1.80, 1.90, 1.7b, 1.55, 1.355, 1,10, 1.045, 

2 .933. . 447. .77. .705. .545. ,bO, ,5b, 1.58, 1.4b, 1.35, 1.445, 1.43, 1.32, 

31. i, 08, , 97, .07, .79,. 72, .bb,.bl,, 57, 1.2b, 1.25, 1,235,1,22,1.20, 
41# 4l,l. 145,1, Ub#a995#«91#,03##7 5b ##b93, #t>4#. 59, 1.05,1. 05,1. 055, 

51  .  j  55 , 1 , 055 # 1 • 07  # 1 . 05 , 1 . 01 5 # • 97  # .925 # • bb  # . 00 , . 74  # . bo , . 03 , . 90 # .90 # 
6, 9  <*  4, , *03, , 905, .9U  7, .91, .91, .90, ,88, ,8b, ,03, .77, .71, .b7, .79, .792, 
7.Ow#.o‘>#,0O4#.tfl#.013#.01b#.02#.022#.0l7#.790#.7o#.74#.o9#l,57# 

01 .79,2. 0,1. 0,1, 52# 1.29,1. 12#. 905#. 80# .  795#, 72# .bb5#.b2# .575# .535# 
91. 5', l,b0#2,0#1.72#1.4tt#l, 29,1, 13,1.0#. 9O#.01#.74#.o0,.o3#.50i# 
1,543,1. 57.1. 55,1. 52,1. 40,1.44,1. 29#l.lb#l.O3#.93#.035#.7b#./O# 

2. o43,.o0, .5b, 1.27,1.30, 1.30,1.20, 1.24, 1.19, 1.145,1.05, ,9b, ,07, .00, 
3.735#, ob#.o25#.50#l. Ob #1,055,1. 05,1.04,1. 02#. 995#.90#, 97#. 97#.90# 
4.03# . 7 o 5 # • 7 u # .055# « b 1 # .905# .90 # .005# .0  75  # . U b 5 # ,0b# .050  # . 055#  .053  # 
5.O5#.O4#.705#.735#,b9#.b5#.795#.703#.70#.77#,7b2#.7o#.757#.75#.75# 
b  •  7 5 #  .  758  #  «  7 b  #  ,  7 7  #  ,  7 4 #  .  7 2/ 

UA1A  bAN4  /  .250# .350, .450# .4b5# .475# .403# .40b# .490, .30b, 

1 • 4  7  4  # , 472 • • 47  2  # . 400 , , 487  # . 490  # . 495  # . 5  # , 5  # .5  # .5  # .5  # ,5 , . 5ol # , 50 1 # 

2. 505..  54 ..54 #.54 ,.54,. 54 ,.525 #,520 ,.505 ,.54  2 ,.500 ,.500,. 580 ,.500, 
3.39u,<558,.509,,548,.5o5,.b25,,b25,,b25,.b25,.b25#.515,.555,.595, 
4,b33, .b73, .b73, ,b73, ,b73, .145, .295, .450, .*75, .405, .490, .495, .495, 
5. 2b /,. 350,. 43 7,. 47 3,. 483,. 487,. 492,. 493,. J80,. 445,. 445,. 445,. 480, 
b. 502, .50 2, .502, .395, .440, .401, .481, .481, .532, .530,. 523, .4U5,. 433, 

7. 470..  525.. 525.. 525.. 5 70..5b0,. 41 3,. 457,. 503, .550,. 570, .570, .570, 
0.bl5,.425,.45O,.477,.5lO,.5bO,.b25,.b25,.b25,.o7O,.2bO,.4bO,.47 7, 

9.445. . 490# .493, .493, .1 7  7 , .27 b , . 3 b 5 , .457 , .475, .402, .405, .405, . 2 0 5 , 
1. 3ow,. 30b,. 41 5, .480,. 405,. 405,. 48 5,. 305,. 355, .400,. 433,. «o5,. 511, 

2. 510. .  5U5, . 317, .323,. 32b, .4b0,. 407, ,501, .543, .537, .330, .320, ,310, 

3. 425..  505.. 523.. 550.. 505.. 33 7.. 315.. 300.. 423.. 50b,. 5b0,. 500,. b05, 
40.00, . 330 , .445, . 4b  5, ,475, .403, .405, ,485, .005, .310, ,4o0, .440, .455, 
5. 4b5, .470, .470, ,lb5, .290, .380, .425, .448, .458, ,4b/, .475, .193, .300, 

6. 390. .  433. .457. .470. .4/7, , 477, .200, .33u, .423, .4b 5, .4v5,, 493, .49b, 

7. 495..  200.. 330.. 423.. 4b5,. 485,. 493,. 49b,. 495,. 200,. 330,. 423,. 4b5, 

8.405. . 493. .49b, .495/ 

XaAN»SW0T( ABSC 1 .-XN4*2> 7 
7aAR*TAN (S*bEP/57 .295781 
l*  Th 

IF ( A.GT , 7 . )  WHITE  ( b , 900 )  X 
lFU.CT.l.)  WRITE  (b,9u2>  Z 
IFlXW.Ct.1.0)  CO  TO  1 

CALL  TABLE3  ( 2 , CLAR , BAR 1 , X , XII , 15 , 3 , 7 , 7 T 1 , 7 , 3 , Z , ZT1 , 4 , 3 ) 

IFU.GT.4.J  WRITE  (6,9011  V 

CAM,  1ABLE3  ( 2,XCPCB,BAR2,X, *12,8, 3 , 7, KT2,5, 3, Z,ZT2, 4, 31 

CLAaCLAK*AR 

RETURN 

1  CALL  IA8LE3  (2, CLAK, BAK3, X, XII, 15, 3, Y,VT1,7, 3, Z,Zn, 4, 31 
IFd.CT.b.J  WRITE  (b, 901 1  7 

CALL  TAULE3  (2 ,XCPC8,BAR4,X,XT2,8, 3, 7 ,VT4, 7, 3,2,112,4,31 

CLA«CLAR*AK 

RETURN 

990  FORMAT  (IX, 'EXTRAPOLATION  REQUIRED,  1 10HAR • BETA  a  F10.4) 

901  FORMAT  ( IX, • EXTRAPOLATION  REQUIRED,  ' 1 1HAR ' TANSw  a  F10.41 

902  FURMAT  ( IX,  •  EXTRAPOLATION  REQUIRED,  TAPER  RATIO  a  'MO. 41 
END 

AUiROUTlNE  MCDEVT  (FMACN, AR ,CLA) 

FIN  NORMAL  FORCE  FROM  NACA  RfcPT.  1253  (MCDEVITT1 
RECTANGULAR  ONLI 

DIMENSION  TFMACH( bl , TAR ( 7 1 , TCLA( 7, b, 1 ) 

DATA  (TFNACH(I1,  lai,6)  /  .  0 ,  .  5  ,  .  •  ,  .  9 , 1 .  ,  1 .  I  / 
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0A1A  (tARH?,  I»1  ,7)  /  .U , . 5, 1 . , 1 .5 , 2 . ,  J.TiT 
DlfA  TCLA  /.O,.Ol4,.U29,„04,.u49,.0b3,.U75, 

♦  •0,.014,«029».04,.049,,Ub3,.075,.0,.U14,.U29,.U4l,*US2,.OB),.U/S 

♦  •  0 , tUlbf «  0  3 , .043# tOSbf iObl , •  1 U  7  , , 0  , *018*  »014| . 0  5  , #Ub4i . 0  8  b  , *099# 

♦  .0, .018, .034# .048, ,0b, .078, .089/ 

call  IA»lE2  (2#CLA,TCbA,AK,TAR,7,2,FMACH,TFMACH,b,2) 

CLAaCL A<S7 . 3 

Rfclunn 

END 

SUenuUUNt  FAC  I  OK  (bU*.Rb*,CK»#AHX,XLN,DELT*,bOT,RbI, IRT,ALI, ACPI, 
C  ALPHA, XL, F*S*) 

DIMENSION  f 01 4 i 

DETERMINATION  Ob  TAIL  INTERFERENCE  FACTOR  U)  bY  STRIP  I’rtEURY 
AND  oLEnDER  BUOY  THEORY  FROM  NACA  lil)7  -  APPENDIX  b 
USE  I R  IS  IN  E0UAI10N  (S5)  10  CALCULATE. 

LIFT  ON  1 AlO  SECTION  DUE  TU  *IN0  VORTICES 

*1*G  *  UK'lb  A  LATERAL  POSITION  IS  SET  AT  .  7854  1SEE  EJUATIJN  J9) 

8u«  x  EXPOSED  SEMI-SPAN  OF  *1nG 
KB*  «  KAU1US  OF  bOD X  AT  *If»G 
CRN  a  ROOT  CHORD  OF  KING 

AH*  x  DISTANCE  FROM  *1nG  LEADING  EDGE  TO  *ING  HINGE  LiNE 

XL*  X  DISTANCE  FROM  NOSE  OF  SUUX  TO  *ING  LEADING  EDGE 

i  DELT »  x  «1nG  DELFECT10N  ANGLE  (SEE  PAGE  11)  -  DeG 
BUI  x  EXPOSED  SEMI-SPAN  OF  TAIL 
RBI  X  RADIUS  OF  BODY  AT  TAIL 
IRI  X  TAPER  ratio  OF  TAIL  1CT/CR) 

ALT  x  DISTANCE  FROM  NOSE  OF  BODY  TO  TAIL  LEADING  EDGE 

ACPI  x  CENTER  OF  PRESSURE  OF  TAIL  (MEASURED  FROM  LEAD  I EDGE  UF 

ROOT  CHORD) 

ALPHA  x  ANGLE  OF  ATTACK  •  DEG 

XI  x  TAIL  INTERFERENCE  FACTOR 

F *S»  x  RATIO  UF  yQRTEA  LATERAL  POSITION  TO  *1NG  EXPOSED  ot-U-SPAN 
(MEASURED  FROM  BODY  SURFACE) 

F*  s  SPANNISE  LOCATION  OF  «1NG  VORTEX  (FROM  BODY  CE>«  TEkLINE) 

FT  x  SPANXISE  L0CAT1UN  OF  «1NG  VORTEX  AT  TAIL  FIN  STaDun 

hT  X  HEIGHT  OF  XING  VORTEX  ABOVE  TAIL  FIN  STATION 

HI  x  EOLATION  41,  PAGE  li 

RAOXS7.295/B 

Pl«3. 14159 

HI»-  (CRN-XHn  )  *S1N  (  0£LT*/K AD  j  4(XLT4XCP1-AL*-CR*)*SIn(ALPi1A/RAu) 
ROSxxRBn/ (HBw+BON ) 

C1X1.-R0SM 

C2«1.4RQSb**2 

iC3xl.-ROS***2 

FnS*x((P1/4.) *C3-R0S*4(C2**2/ (2.*C3))*ASIN(C3/C2))/(2.*C1) 

.  FNxf nSn*B0n  ♦  Rbx 

FTxr«-RbM4HbT 

8«BUT4RbT 

RSxRbT/S 

HSxhI/S 

FSxFT/S 

1F((FS.E0.1.).AND.(HS.EO.O.))  FSX.99 
rS2x(F5*RS**2)/iFS**2  ♦  HS**2) 

HSI*(HS*RS**2 )/ (F5**2  4  HS**2) 

1X1 

•0  TL( 1 )>( 1 ( 1 .0-RS*TKT )-FS* ( 1 .O-TRT ))/(2.0*( 1 . U-RS) ) ) 
TL(I)xTLlJ) *ALOG( ( HS**2  4  IFS-l ,0)**2)/(HS**2  4  (FS-RS)**2)) 
TLll)XTL(l)-ll.-TRT) 

Xf(HS.EO.O.O)  GO  TO  «1 

TL(I)*TL(l)-(( 1 .-TRT) /( 1 .-RS) )*HS* ( ATAN( (F5-l.)/riS) 

1  -ATARI IFS-R5)/HS) ) 

•1  1»1  4  1 

60  TO  (12, 82, 83, 84, 85), 1 
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*4  FS9-m 
GO  XL  «u 
91  FSaFSl 
HS*nSl 
GO  To  90 

84  FSa-FSl 
GO  JO  HO 

85  Xia(2.0/(1.0  ♦  THXJ  )*(lL(l)-TL(2)-'fL(l)*TL(-*>) 

RElURh 

END 

SUbHOUlINE  FNCKU5  ( FMACH ,  AK ,C1  ) 

THIS  SUBROUTINE  CALCULATES  THE  CUEFF.  OF  NUN-LINEAR  CunTKIBUilUN 
TO  LlKT,  Cl. 

TAKEN  FROM  A  RARER  PRESENTED  AT  ThE  SIMM  U.  o .  "tA  VAL  S(*POSuIM 
ON  AERQBALL1ST1C5,  I  NUV  19bl. 

AUThOk,  a.  f.  lakpros 

note:  THIS  THEURY  IS  VAulD  only  FOR  SUHS0N1C  F LJw • 

(Cl)  IS  REDUCED  UUADKAT1CALLY  *  IT H  MACH  NUABEK 
BEYOND  AACH  1.0.  THIS  IS  AN  ARBITRARY  FIX  1U  SMOOTH 
THE  TRANSITION  INTO  SUPERSONIC  FLU*. 

DIMENSION  TC1U0),  TAR1U0) 

DAI A  (lCl(l), 1«1, 101/9.8, 8. 95,6.bb,5. 1,4.45 ,1.75. 1.1, 2. 55, 4,, 1.55/ 
DATA  ITARlll),  1*1,10)  /.O, .5, I., 1.5, 2., 2 .5, 3., 4., 5., b./ 

CALL  1NTERP  C2,1,C1,TC1,AK,TARI,10,3,w1N,maX) 

Cl*Cl*.0UOi 

IF  (FAACH.GT.1.)  C1*C1/FAACH**4 

RE1URN 

END 

I  SUbkuLTlNE  VORTEX  ( BOW , KB* ,CRw , XH* , XL* .DELTA , Phi  * , A* AR , BU t , Rbl , 
j  1  IR1, XLT, ACPI, PHI T, ART, AdPHA.CLa* ,CLAT , BK*B , SK *B , Cb I V , F  *S*) 

!  THIS  SUBROUTINE  CALC.  LIFT  ON  TAIL  SECTION  DUE  TO  *1*0  VORIlCbS. 
j  iHE  METHOD  OUTLINED  IN  NACA  110)  IS  USED. 

ApPbNDlX  b  15  USED  FOR  THE  TAIL  INTERFERENCE  FACTOR  (1). 
INDEPENDENT  ROLL  ANGLES  OF  WING  AND  TAIL  ARE  ALLUnED. 

THE  SUBROUTINE  ACCOUNTS  FOR  ROLL  ANGLES  BX  ADJUSTING  THE  GEOMETRY 
AND  TAKING  COMPONENTS  IN  THE  PITCH  PLANE. 

THE  FINAL  ANSWER  (CLTV)  IS  VECTORED  IN  THE  ANGLE  OF  ATTACK  PLANE 

(OK  PITCH  PLANE)  REGARDLESS  OF  WiNG  OR  TAIL  ROLL  ANGLES. 

BOW  *  EXPOSED  SEMI-SPAN  OF  WlNG 
RBW  a  RADIUS  OF  BOOT  AT  wlNG 
CRN  *  ROOT  CHURO  OF  wlNG 

XHw  •  DISTANCE  FROM  WlNG  LEADING  EDGE  TO  WING  HINGE  LINE 

XL*  *  DISTANCE  FROM  NOSE  OF  BOOT  TO  WlNG  LEADING  EDGE 

DELTw  a  wlNG  DCLFECT10N  ANGLE  (SEE  PAGE  11)  -  DEG 
PHIw  a  ROLL  ANGLE  OF  wlNG  FROM  VERTICAL  (DBG), 

ANAR  a  RATIO  OF  wlNG  AREA ( 1*0  wINGS)  TO  REFERENCE  AREA. 

BOT  a  EXPOSED  SEMI-SPAN  OF  TAIL 
RBI  a  RADIUS  OF  BODY  AT  IAIL 
TRT  a  TAPER  RATIO  OF  TAIL  (CT/CK) 

XLI  a  DISTANCE  FROM  NOSE  OF  BOUY  TO  TAIL  LEADING  EDGE 
XCPT  a  CENTER  Or  PRESSURE  OF  TAIL  (MEASURED  FROM  LEADING  EDGE  UF 
HOOT  CHORD) 

PHIT  a  ROLL  ANGLE  OF  TAIL  FROM  VERTICAL  (DEG). 

ART  a  ASPECT  RATIO  OF  TAIL  (TWO  TAILS  JOINED  WITHOUT  BODY). 

ALPHA  a  ANGLE  OF  ATTACK  •  DEG 

CLAW  ■  LIFT  COEFFICIENT  SLOPE  OF  WlNG,  BASED  ON  wlNG  AREA  (1/RAD) 

CLAT  a  LIFT  COEFFICIENT  SLOPE  OF  TAIL,  BASED  ON  TAIL  AREA  (1/RAD) 

BKWB  a  UPWASH  INTERFERENCE  FACTOR  FOR  wlNG,  DUE  TO  ALPHA. 

SMB  a  UPWASH  INTERFERENCE  FACTOR  FOR  wlNG,  DUE  TO  DELI*. 

SEE  NACA  WO? ,  PAGE  49,  FOR  BKwB  AND  SKWB. 

CLTV  a  LIFT  COEFFICIENT  ON  TAIL  SECTION  DUE  TO  RING  VORTICES, 
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OH^tU  ON  KfcFEKbNCt  AREA. 

htJuHw  *■  ilh  CCtlVJ  6UHMLU  t  OK  ALL  CA5tS. 

Mr  FtCTlNG  EACH  1 AIL) . 

XI  *  t  AIL  INTERFERENCE  ft  ACTOR.  (UNE  *lNG  ACTING  On  jNfc  lAlU) 

H»5*  *  H  A  U  0  OF  VORTEX  LATfcKAL  PUSIIIUN  TU  »1NG  EXPOSED  SEMI-SPAN 

tMfcASUKtO  fr  HON  BOO*  SURFACE) 

OlNENSiUN  PH«lNbl4>,PHlAlLl4),HI(4),Fll4)  ,PH1V(4),RV14)  ,  IL(2) 

UNNbNl,. .  NO.  OF  DING  VuKTEA  AND  TAILS 

Pl«J.14li>92b5JbU 
RAO**/. 2957/9512* 
lF(Pnl*)  7,1,7 
7  lF(PHiw-9U. )  2,  1,2 

1  NVONTS2 

PHNlNL(l)*Fi/X. 

PM*ING(2)b3.*P1/2. 

GO  10  3 
.  2  Ntf  UN  1*4 

•  PNNli'.OU  )«PH1m/HAD 
(  PHn1nG(2)s(PhIn*90.)/KA0 
{  PHNlNGli)*iPHlD4ltiO.)/NAU 
PH»1NG(4)b(PN1ii-*270.)/RAO 
!  3  lF(PHlT)  t),4,tf 

I  B  iF(PHlT*9U.)  5,4,5 

!  4  NlAIL*l 

PHXA1L(1)bP1/2. 

GO  10  b 

5  NTA1L*2 
PHTA1LI 1 )*PN1T/NAD 
PHTAlL(2)*CPHIT+90.)/RAD 

6  CONTINUE 

OHNENT.,..., . . . .  wlNG  VOhfbA  PUiiUUN 

ROSw«Kbw/(RBwH*Uw) 

0*1. -RuS* 

C2«l.fNU5»**2 

jC3*l.*RUSN**2 

F*Sw*(CP1/4.)*CJ-K05r*(C2A»2/(2.*C3))*AS1N(C3/C2))/U.*C1> 
FN«FHS«*BON-»KbR 
FNaf  m-RBn+RBT 
SaBOl+RbT 
KSbNbT/S 
,00  10  K*1,NV0KT 
SlPH««SlNlPHnING(K)  ) 

!C0PMm«CUS(PHn1NG(K)) 

;SIOL*»SlN(OfcLTw/RAO) 

SI ALP*S1N (ALPHA/ HAD) 

HT(H)«-(CRw-XH*)*SlDb*«*AttS(SlPHl») 

C  ♦F«*C0PHNMAbT4XCPT«AbN*CRM)*SlALP 

>FT(K)arN*61PH«4(CRN*XHN)*S10bM*CUPHli*SlPHN/ABS(SlPH«l) 
PH1*ATAN(HT(K)/FT(K)  ) 

PHl»(R)«Pl/2.-PHl 

IFCFl(K) )  9,9,10 

9  PH1V(K)«I,*P1/2.*PH1 

10  RV(K)a3QftT(FT(K)**2  ♦  HT(K)**2) 

ONRCNT . . . .  CALC.  TAIL  iNTERFEK,  FACTUK 

CbTVNO.O 
00  15  L*1 ,2 
00  11  Oal,ATAlb 
00  11  K*l , ft VORT 

fVNV (K)*COS(PHTAIL( J)-PH1 V (A) ) 

H«RV(K)«SlN(PHTAlb(0)-PHlV(X7) 

riar/s 

P«la(rS»RI**2)/(FS9*2  ♦  HS*«2) 
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US**2  ♦  HS**2 ) 

00  12  lsl,2 

TL(l)s(((l.0-RS4THT)-FS*(1.0-tRT})/(2.0*(1.0-RS))) 
TL(1)*IL(X)*ALUG(  (HS**2  ♦  (FS-1.0)**2)/(HS*»2  ♦  (FS«KS)**2)) 

TLU )«tL(l)-t(l.O-rRl)/Cl.0-HS) )♦( (1,0-RS)  ♦  HS*(AIAN((FS-1.0)/HS) 
C  -AiA.n  (FB-K5)/HS)  )  ) 

FSafSX 

12  H5*MSi 
XL  sTLll) 

XLIM  mTLU) 

XIs(2./(1.4TRT))*(XL-AL1M)*SXN(PHN1NG(K)  )/AbS (SI N ( PH*1 NG{ K )  )) 

CiNMENl. .................. ....... .  CALCULATE  CbTv 

GAMOSVa (BKNB* ( ALPHA/RAU  ) * ABS ( SI N ( PHN 1 NG ( A ) ) ) tSK * B*UEL T*/K AO  ) 

C  *CLA«*A*AR*AbS(S1N(PHm1NG(K)  )  )/(4.*FnSn*BCJ*  ) 

CLPAHT*Xl*CLAX*AbSlSlN(PrtlAlLU))  )  *»OX*iiAAJoF/(Pi*AK,i  ) 

CLTvaCLTV+eLPAKl 

11  CONTINUE 

00  li  J*1 , NTAlb 

li  PHTAXL(J)*PI-PHTAXL(J) 

DO  14  K*1 , NVOHI 
PHnING(K)b2.*PI-PHi»1NG(K) 

14  PHlV(K)a2.*Pl-PHlV(K) 

15  CONTINUE 
RETURN 

END 

SUBROUTINE  INTSLB  ( N , F* ACh , BETA ,  DM , RB , CR , TR , BU , AH ,CbA , 

1  bKNB.bKBM.SKNB.SKBM) 

nXNG-BOU*  A NO  TAIL-BUD*  INlEKt ERENCE  FACTORS  FROM  CHARTS  1,4,5 
NACA  1307 

N  »1  At  TERB0D1 

N  s  2  NO  AFTERBOU* 

FA ACH  *  FREE  STREAM  MACH  NUMBER 
BETA  s  SUBSONIC  ON  SUPERSONIC  BETA 

BN  s  BETA  *  COTANGENT  uF  LEADING-EDGE  SmEEP  AmGLE 

RB  s  RADIUS  OF  BUD* 

CR  a  ROOT  CHORD 

TR  a  TAPER  RATIO  (CI/CR) 

BO  a  EXPOSED  SEMI-SPAN 

AR  a  EXPOSED  ASPECT  HAI10 

CLA  a  LIFT  COEFFICIENT  SLOPE  (PER  RADIAN)  -  BASED  ON  NXNG  OR 
TAIL  EXPOSED  AREA  (2  PANELS) 

BKNB  a  BIG  (K)  NING-BQD*  INTERFERENCE  (CHART  1) 

BKB*  a  Bltt  (K)  BOOT-MXNG  INTERFERENCE  (CHARTS  1.  4.  5) 

SKNB  a  SHALL  (K)  »1nG-B0DT  INTER! EKENCE  (CHART  X) 

SKBN  a  SMALL  U)  BOO*-«XNG  INTERFERENCE  (CHART  X) 

DIMENSION  TBKnB(  1 1 )  . TBKBn  ( 1 1 )  , TSKNB ( 1  1 )  .  TSKBN (  1 1  )  , TKOS(  1 1 )  . 

1  TBKC(7J,TBRCN(7).fBNUV).TBNN(«).TKBNC(7.19).rKBNCN(7.B) 

DATA  (TBKNB(l)  .  1st. 11)/ 

X  i. 0, 1. 08, 1. 165, 1. 255, 1.352,1. 4b, l.SSS, 1.663, 1.7 72, 1.8*5. 2.0/ 

DAtA  (TBKBM(l).  lal.ll)/ 

X  0.0, 0.1 3, 0.28, 0.44, 0.61,0. BO, 1.0, 1.22, 1.4b, 1.70,2.0/ 

DATA  (TSKMB(X),  lal.ll)/ 

X  1.0, .963, .943, .93$, .935, .939, .946, .958, .970, .90S, 1.0/ 

DATA  (TSKBN ( 1 ) ,  lal.ll)/ 

X  0.0, 0.11,0, 21, 0,31b, 0.41b, 0.515, 0.61,0. 70S, 0.002, 0,40, 1.0/ 

DATA  (TROSU),  lal.ll)/ 

X  0.0, 0.1, 0,2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, 1.0/ 

DATA  (TBRC(i),  1*1,7)/  ,1, .6, 1.0, 1.4. 2. 0,2. 8, 4.0/ 

DAtA  (TBRCN(l),  Xal,7)/  ,2* .6. 1 .0, 1 .4,2.2, 2.8 , 4.0/ 

DATA  (TBM(l),  1*1 .18)/  , 1 , ,2 , , 3, • 4, *S , ,6 , , 7 , , 8 , .9 , 

X  1 .0, 1 .2, 1 .6, 2.0. 3.0, 4.0, 7.0, 10 ,0,99 .9999EF 30/ 

DAtA  (TBMNl 1 ) ,  lsl.l)/  .2, ,4, .6, .8, 1 ,0, 2 .0,4 .0,99. 9999C+ 30/ 

OAtA  TKBNC/ 
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X 

.0331,  .349b* 

.2665* 

.1832,  .1333# 

X 

1.2b#  . b 1 6  3  # 

.6664* 

jfwrgpgTFPjHl 

,4332,  ,3998* 

X 

1.806,  1.2b# 

1.05*  . 

933*  .8163*  .6831#  .8331* 

X 

<•332*  1.060 

*  1.399# 

1.25*  1.0U99* 

.9496*  .8497* 

X 

2.849*  2.010 

#  1.749, 

1.549*  1.366* 

1.199*  l.yb, 

X 

3.249*  2.316 

*  2.016, 

1 • 8l6  *  1.633# 

1 , 460  *  1 . 2  3  3  # 

X 

3.61b*  2.632 

#  2 . 266  # 

2.066*  1.893* 

1.633*  1.433* 

X 

3.898,  2,862 

#  2.499, 

2,266,  2.0*  1 

.799,  l.bbb. 

X 

4.21b*  3.149 

*  2.716* 

2.466*  2.199* 

1.9oo,  1.733* 

X 

4.41b*  3.31b 

*  2.882* 

2,599*  2.332# 

2.083*  1.832* 

X 

4.781*  3.61b 

*  3.149* 

2.866*  2.51b* 

2.249,  2.010* 

X 

6.181*  3.998 

#  3.515* 

3.21b*  2.849* 

2.516,  2.232* 

X 

b.blb#  4.24b 

.  3.749* 

3.432*  3.U82* 

2.732*  2.382* 

X 

6.914*  4. bob 

*  4.132* 

3.798,  3.399* 

3.049*  2.682* 

X 

6.214*  4.91b 

,  4.382# 

4.ulb,  3.632* 

3.232*  2.866* 

X 

0.647*  b.26b 

*  4.681* 

4.31b,  3.898, 

3.499,  3 . 0»2 

X 

6.864*  b.431 

*  4.848, 

4 ,46b  #  4.048, 

3.032,  3.21b* 

X 

7.23*  b.848# 

5.231, 

4.831*  4.365* 

3.932,  3.482/ 

DATA  TKBWCN/ 

X 

.  9bb3  #  .5442 

*  .3685. 

.2632*  .1755, 

.1229*  .0877* 

X 

1.86*  1.14* 

•  7  3  7 ,  b  333  p  *24b 

#  •  1  7  6  * 

X 

2.b97 #  1.509 

#  1.017, 

.737, .456*  ,351*  .281* 

X 

3.123*  1.772 

*  1.176, 

.84,  .544,  .439*  .316* 

X 

3.65*  2.264* 

1.457* 

1.053#  .684,  . 

544  #  .368* 

X 

4.56#  2.861# 

1.842* 

1.316*  .842*  . 

666*  .491, 

X 

5.247,  3.287 

#  2.088* 

1.527*  l.UOO, 

.789,  .544, 

X 

6.107*  3.879 

*  2.51* 

1.772*  1.193* 

,912*  .649/ 

KOS*Kb / ( Rb  +  BO) 

IFtFMACH.LE.l.OJ  GO  TO  1 
IF(BM.GT.99,9999E*30)  8M*99 . 99996+ 30 
CON*bfclA*AR*ll.O  ♦  TR)*(1.0  ♦  1.0/  BH) 
if (CUN. 0£. 4.0)  GO  TO  1 
BKC*2  » 0*BETA*R0/CR 
IF  (  6KC . GT .4.0)  BRC*4.0 

C0N2abeiA*CLA4(TR  ♦  l ,0 ) * ( l .0/K05  -  1.0) 
lF(N.GT.l)  GO  TO  3 

CALL  TABLE2  l 2 , ANR , TKBwC , BRC , 1BRC , 7 , 3 , BM , Tb* , 1 b , 3 ) 

BKBWSANR/C0N2 
GO  TO  4 

3  CALL  TABLE2  ( 2 , ANR , TKBwCN , BRC , TBRCN , 7 , 3 , BM , TBMN , 8 , i ) 
BKbP*ANR/C0N2 

GO  TO  4 

1  CALL  INTERP  (2,l,BKBw,TSK8N#ROS,TROS,ll,3,MlN,MAX) 

4  CALL  INTERP  ( 2 , 1 , BKwB , TBKWB, R05 , TRU5 , 1 1 , 3, MN , RAX ) 

CALL  INTERP  12*1* SKBW* ISKB*  *RQS#TROS ,ll#3,HIN,MAX) 

CALL  INTERP  (2,i,SKWB,TSKWB,R0S,TR0S, 11*3, MIN, MAX) 

RETURN 

ENO 

SUBROUTINE  INTERP  UXTRP , LMT , t , YT , X , XT , NX , NPA # MINX , NAXX ) 

TABLE  LOOK-UP  ROUTINE  FOR  1  INDEPENDENT  VARIABLE 
WITH  SIMPLE  INTERPOLATION  UR  EXTRAPOLATION. 

IXTRP  ■  1  EXTRAPOLATE  IF  NECESSARY, 

IXtRP  ■  2  EXTRAPOLATION  NOT  ALLONED,  TAKE  THE  LASItQR  FIRST)  POINT 
LNT  «  1  PROGRAM  MILL  DETERMINE  SUBSCRIPT  RANGE  (MINX  TO  MAXX). 

LMT  «  2  PROGRAM  ASSUHS  THAT  MINX  AND  MAXX  SUBSCRIPTS  ARE  KNORN, 

X  ANSWER  UK  OEPCNDENT  VARIABLE  CORRESPONDING  TO  INPUT  (X). 

YT  TABLE  OF  DEPENDENT  VARIABLES  CORRESPONDING  TO  XT. 

X  ARGUMENT  OR  INDEPENDENT  VARIABLE  FOR  WHICH  ANSWER  (Y)  WILL 

BE  DETERMINED. 

TABLE  OF  INDEPENDENT  VARIABLES. 

NUMBER  OF  POINTS  IN  XT* YT 


XT 

NX 


MPA  NUMBER  Up  POINTS  10  Dfc  USED  FOB  1  i*  )'  EHRUL  A XiUM  . 

MIMA  MINIMUM  XI  SUBSCRIPT  USED  FOR  INTERPOLATION. 

MAXX  MAXIMUM  XT  SUBSCRIPT  USED  FOR  INTERPOLATION . 

DIMENSION  XT (NX J , X  I ( NX J 
IF  UXTRP.LU.1J  GO  TO  140 
faiTU) 

IFU.Lt.XTU J J  RETURN 
X*U(NX) 

1F(X.gE.XT(NX)  J  RETURN 

140  IF(LMX.EU.2)  GO  TO  1J0 

CALL  LIMIT  IX, Xi  ,NA, NPX, MINX, MAXX) 

110  ysTKMlMXJ 
IF(MINX.EU.MAAX)  RETURN 
**0. 

UU  120  JaM 1NX  > MAXX 
P*I. 

DO  UU  laM 1NX ,  MAX  A 
j  1FU.E0.JJ  GO  TO  110 

i  p*p*ix-xruj j/uuj)-xtujj 

I  UU  CONTINUE 

I  1 20  Tax  ♦  IIU)f? 

i  RETURN 
END 

SUBROUTINE  TABULA  ( IXTHP , w , WT , X , XT , NX , NpX , X , XT , NX , NPX) 

TABLE  LOOK-UP  ROUTINE  FOR  2  INDEPENDENT  VARIABLES. 

1XTRP  a  1  EXTRAPOLATE  IT  NECESSAKX . 

1XTRP  a  2  EXTRAPOLATION  NOT  ALLOWED.  TAKE  THE  LAST ( OR  FIRST)  POINT, 
j  w  a  ANSaER, (DEPENDENT  VARIABLE  CORRESPONDING  TO  INPUTS  XU) 

*T  a  TABLE  UF  DEPENDENT  V AR1BALE  CORRESPONDING  TO  XT, XT 

ntu.j)  Increment  subscripts  left  to  right  when  loading 

X  a  the  argument  UR  INDEPENDENT  VARIABLE  X 

XI  a  TABLE  Ut  1N0EP «  X  VALUES  (MUSI  BE  IN  INCREASING  ORDER) 

NX  a  NUMBER  UF  POINTS  IN  XT 

NPX  a  NUMBER  U|  POINTS  TO  USE  FOR  X  INTERPOLATION 

X  a  THE  ARGUMENT  OR  INDEPENDENT  VARIABLE  X 

XT  a  TABLE  UF  INDEP.  X  VALUES  (MUST  BE  IN  INCREASING  ORDER) 

NX  a  NUMBER  OF  POINTS  IN  XT 

NPX  a  NUMBER  OF  POINTS  TO  USE  FOR  X  INTERPOLATION 

DIMENSION  XT(NX),XI(NX),WT(NX,NX),  B(20) 

CALL  LIMIT  (X, Xf , MX, NPX, MINI, MAXX) 

CALL  LIMIT  (X, XI, NX, NPX, MINX, MAXX) 

1  DO  42  JaMlNX , MAXX 

CALL  INTERP  (IXTRP,2,B(0) ,WT(1,J) , X, XT, NX, NPX, MINX, MAXX J 
42  CONTINUE 

CALL  INTERP  ( IXIRP, 2, W ,B, X , XT ,NX ,NPY , MINX ,MAXX ) 

RETURN 

END 

SUBROUTINE  TABLEj  ( IXTRP, w, wT,X,XT,NX,NPX, X , XT , N t , NPX , Z,2T , NZ.NPI) 

TABLE  LOOK-UP  ROUTINE  FOR  J  INDEPENDENT  VARIABLES. 

IXTRP  a  I  EXTRAPOLATE  IF  NECESSARX. 

IXTRP  a  2  EXTRAPOLATION  NOT  ALLOMED.  TAKE  THE  LAST(OR  FIRST)  POINT. 
■  a  ANSWER, (DEPENDENT  VARIABLE  CORRESPONDING  TO  INPUTS  X,X,I) 

NT  a  TABLE  OF  DEPENDENT  VARIBALE  CORRESPONDING  TO  XT,X?,ZT 
wT( 1 , J ,K)  INCREMENT  'UBSCR1PTS  LEFT  TO  RIGHT  WHEN  LOADING 
X  a  THE  ARGW  ENT  OR  .DEPENDENT  VARIABLE  X 

XT  a  TABLE  Or  .a  EP  X  VALUES  (NUST  BE  IN  INCREASING  ORDER) 

NX  a  NUMBER  OF  : DINTS  IN  XT 

NPX  a  NUMBER  OF  POINTS  TO  USE  FOR  X  INTERPOLATION 

T  a  THE  ARGUMENT  OR  INDEPENDENT  VARIABLE  X 
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II  s  iAbLE  0F  INOEP.  1  V  AbUtS  (MUST  bb  IN  InCKEASI  -Hi  UHUEh) 

M  s  NUMbER  Uf  POINTS  IN  IT 

NPX  x  NUMbER  UF  POINTS  1U  USE  FUR  K  INTERPOLATION 

Z  x  Iht  AKCUMbNl  UN  INDEPENDENT  VARIABLE  Z 

ZT  *  TABLE  UF  INDEP.  L  VALUES  (MUSI  BE  IN  INCREASING  UNDER ) 

Si  s  NUMbER  UF  POINTS  In  U 

N PZ  *  NUMbER  UF  POINTS  Io  USE  FUN  Z  INTERPOLATION 

DIMENSION  XT(NX),  *T(NK)#21(NZ),NT(NXfNX,NZ>,  b(2U),A(10) 

CALL  LIMil  (Z,ZI,NZ,NPZ,MINZ,MAXZ) 

CALL  LIMIT  ( X , 11 , Hi  ,  N  P  X f  M I N I ,  M  A  X  X ) 

CALL  LIMIT  (X,AT,NA,NPA,M1NA,MAAX) 

DO  41  AXM1nZ,MAAZ 
DO  42  Jsm I NX  r  M A A I 

CALL  InTEKP  (IAl'RP, 2, b(U), *1(1, U,K),A, XI, NX, NPX, MINA, MAAA) 

4k  CONTINUE 

CALL  IN  TEMP  ( I X I R P , 2 , A ( N ) , b , X , IT , NT , NP X , M 1 n t  ,  M A X I J 
41  CONTINUE 

CALL  InTEKP  (lXTRP,k,*,A,Z,ZT,NZ,NPZ,MINZ,MAA2) 

NElUNN 

END 

SUbNOUT 1NE  LIMIT  (X,XI,NA,NP  ,MINX,MAXX) 

THIS  SubKOullNE  MILL  FIND  THE  MINIMUM  AND  MAXIMUM  SUbSCKIPTS 
ION  NANCE)  TO  bb  CONSIDERED  FOK  INTERPOLATION. 
ulMENSlUN  Al(NA) 

npx*np 

IF(NPX.GT.NX)  NPX*NX 
DO  kb  1*1, NA 
IF(XiU)-X)  2b, 22,21 
kb  CONI INUE 

.  CHEATER  THAN  maA  SUBSCRIPT 

24  MAXX*NX 
M1NX*NX-NPX+1 
RETURN 

. . . . . .  NlTHlN  RANGE 

kl  MlNX*I-NPX/2 

MAXA*MlNX+NPX”i 
1FCMAXX.GT.NX)  CO  TO  24 
If  ClNX.GE.l)  RETURN 

. . . . . .  LESS  THAN  MIN  SUBSCRIPT 

I  MIN  X* 1 
MAXXBNPX 
I  RETURN 

. . .  NO  INTEKP  NECESSARX 

2k  M 1 NX*1 

MAXX*1 

RETURN 

END 

ABEND 
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